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ABSTRACT 

We present the results of a Hubble Space Telescope NICMOS imaging survey of 22 T-type field brown 
dwarfs. Five are resolved as binary systems with angular separations of 0'.'05-0'.'35, and companion- 
ship is established on the basis of component F110W-F170M colors (indicative of CH4 absorption) 
and low probabilities of background contamination. Prior ground-based observations show 2MASS 
1553-f 1532AB to be a common proper motion binary. The properties of these systems - low multiplic- 
ity fraction (llt3% resolved, as corrected for sample selection baises), close projected separations (p 
= 1.8-5.0 AU) and near-unity mass ratios — are consistent with previous results for field brown dwarf 
binaries. Three of the binaries, 2MASS 0518-2828AB, SDSS 0423-0414AB and SDSS 1021-0304AB, 
have components that span the poorly-understood transition between L dwarfs and T dwarfs. Spectral 
decomposition analysis of SDSS 1021-0304AB reveals a peculiar flux reversal in this system, with a 
T5 secondary that is ~30% brighter at 1.05 and 1.27 fjjn than the Tl primary. This system, 2MASS 
0518-2828AB and SDSS 1534-1-1615 all demonstrate that the J-band brightening observed between 
late-type L to mid-type T dwarfs is an intrinsic feature of this spectral transition, albeit less pro- 
nounced than previously surmised. We also find that the resolved binary fraction of L7 to T3.5 dwarfs 
is twice that of other L and T dwarfs, an anomaly that can be explained by a relatively rapid evolution 
of brown dwarfs through the L/T transition, perhaps driven by dynamic (nonequilibrium) depletion 
of photospheric condensates. 

Subject headings: stars: binaries: visual — stars: fundamental parameters — stars: individual (SDSS 
J042348. 57-041403.5, 2MASS J05185995-2828372, SDSS J092615.38-h584720.9, 
SDSS J102109.69-030420.1, 2MASS J15530228-f 1532369) — stars: low mass, 
brown dwarfs 



1. INTRODUCTION 

Multiple star systems are of fundamental importance 
in the study of stellar populations, and by inference much 
of galactic and extragalactic astrophysics. These sys- 
tems remain the predominant outlet for the direct mea- 
surement of individual stellar masses, either through the 
detection of or bital motion or microlensing techniques 
l)An et alJl2002D . Eclipsing binaries also enable measure- 
ment of stellar radii. The properties of and interactions 
between the components of multiple star systems are fun- 
damental to the phenomena of cataclysmic variables, X- 
ray binaries. Type la supernovae, planetary nebulae and 
several classes of peculiar stars. Indeed, the creation of 
multiple systems is inherent in the star formation pro- 
cess itself. Measurement of multiplicity statistics — the 
binary fraction, mass ratio distribution and separation 
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distribution — provide key empirical constraints on star 
formation theory. The formation and character of planet- 
forming debris disks around young stars can be modu- 
lated by the presence of companions. Finally, coeval bi- 
nary systems provide a unique control environment for 
studying the detailed physical properties of individual 
stars, yielding insight on the general characteristics of a 
stellar class. 

Multiple systems have been particularly useful 
in the study of brown dwarfs, stars with insuffi- 
cient mass to sustain core hydrogen fusion ijKumail 
Il962t Havashi & Nakano 1963). Indeed, many of 
the first brown dwarfs to be identi fied are mem- 
bers o f nearbv multip l e svst e ms ( Becklin fe Zuckermgiil 
19881 iNakaiima et al.l 119951 lOppcnhcimer et al. 199E 



Rebolo et al .119981) . Over the past few years, high resolu- 



tion imaging and spectroscopic surveys of very low mass 
(VLM; M < O.IMq) stars and brown dwarfs in the field 
and in young stellar clusters have revealed roughly 75 
binaries (cf. Burgasser et al. 2006b^), with intriguing re- 
sults. The resolved VLM binary fraction (the frequency 
of binary systems in a given sample of stars) is ~10-20%, 
significantly lower than the binary fractions of solar-type 
stellar systems (~65%; e.g. Abt & Levy 1976; Duquen- 
noy & Mayor 1991) and early-type M stars (~30-40%; 
e.g. Fischer & Marcy 1992; Reid & Gizis 1997; Delfosse 

^ A current list of known VLM binaries is maintained 
by N. Siegler at the Very Low Mass Binaries Archive, 
jhttp : //paperclip . as ■ arlzona. 6du/$\ siin$nsiegler/VLM_binaries/ 
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et al. 2004), indicating a decline in th e binary fraction 
with later spectral types l)Bouv et alJ 12006'). The re- 
solved binary fraction is likely a lower limit to the true 
binary fraction due to the existence of unreso lved, closely 
separated systems fMa xted fc Je ffries' 2005V This pos- 
sibility is an important consideration for low-mass sys- 
tems, for while the distribution of separations of F- 
through M-type stellar pairs is quite broad, ranging 
over 0.1 AU to 0.1 pc, >90% of all kn own VLM bina- 
ries ha ve projected separations <20 AU l|Burgasser et al.l 
l2006bj) . with maximum separations scaling wit h total 
system mass l)Close et a l. 2003; Buraasscr ct al. 20033, 
however, see Luhman 2004 and Billeres et al. 2005). 
The mass ratio distribution of resolved VLM binaries 
is also distinc t, peaking sharply at q = M2/M1 
1 (IReid et al.l ^OOlj Bouv et al. 200l IBureasser etall 
l2003cj) . in contrast to the relatively flat mass ratio dis- 
tributions of stellar systems (e.g., Mazeh et al. 1992). 
These properties have led researchers to suggest that 
VLM stars and brown dwarfs may form via a different 
mechanism than stars (e.g., Bate, Bonnell & Bromm 
2002), although this idea remains controversial (e.g., 
Luhman 2004). More concretely, astrometric and spec- 
troscopic foUowup of VLM binaries ha ve provided the 
first brown dwarf mass measurements (Basri & Martin' 
K^99; Lane ct al. 2001; Bouv et al. 2004; Brandncr ct al. 
12004; Zap atero Osorio et al. 2004; Stassun ct al. 2006), 
important empirical constraints for theoretical evolution- 
ary models. 

In order to constrain the binary properties of brown 
dwarfs in greater detail, to identify new systems useful 
for mass measurements, and to search for very low lu- 
minosity brown dwarf companions, we have conducted a 
high resolution imaging survey of 22 T dwarfs using NIC- 
MOS on the Hubble Space Telescope (hereafter, HST). 
T dwarfs are the lowest luminosity (L < 3xl0~^ L©) 
and coldest (Te// < 1400 K; Golimowski et al. 2004) 
brown dwarfs currently known. They are distinguished 
by the presence of strong H2O an d CH4 absorption bands 
in their near infrared spectra ijBurgasser et al.l l2002cl: 
ICeballc ct al. 2002) and the absence of photospheric con- 
densates that dominate warmer L dwarf atmosphere s 
^Marlevet aLII 996HTsuii et al .11 99611 A Hard et aLl200lD. 
We have identified five binaries in our sample, of which 
three have well-resolved components allowing detailed 
characterization of their empirical properties. 

Observations are described in § 2, including the sam- 
ple composition, observing strategy and data reduc- 
tion. In § 3 we present photometric results, including 
color/spectral type, photometric conversion and bolo- 
metric correction relations. In § 4 we describe point 
spread function (PSF) fits to our resolved sources, and 
determine sensitivity limits for putative faint compan- 
ions. Detailed analysis of individual systems is given in 
§ 5. In § 6 we provide an updated assessment of the mul- 
tiplicity properties of field brown dwarfs, including the 
overall binary fraction, separation distribution and mass 
ratio distribution. In § 7 we examine what currently 
known brown dwarf binaries reveal about the poorly un- 
derstood transition between L dwarfs and T dwarfs. Re- 
sults are summarized in § 8. 

2. OBSERVATIONS 



2.1. T Dwarf Targets 

Observations presented here incorporate data from 
two HST programs, GO-9833 and GO-10247 conducted 
during Cycles 12 and 13, respectively. The first pro- 
gram targeted 22 T dwarfs identifi ed in the Sloan 
Digitial Sky Survey iYork et all 1200(1 hereafter SPSS) 
and the Two Micron All Sky Survey ISkrutskie et alJ 
l20f)6l hereafter 2MASS), spanning the full range of 
T spectral types (TO to T8), including the low- 
est luminosit y brown dwarf so far i dentified, 2 MASS 
0415-0935^° feurga,sser et a,].l 12002^ IVrha et a,1.l I200I 
IGolimowski et alJ l2004|) . Twelv e of these source s 
have m easured parallaxes from Dahn ct al. (2002jk 
iTinnev. Burgasser & Kirkoatrick (2003); and Vrba et all 
(Hn^); and 17 have proper motion measurements (al- 
though several did not at the time HST images were 
obtained; see below). Program GO-1 0247 targeted th e 
peculiar Tl dwarf 2MASS 0518-2828 ^Cruz et alJl200l . 
a source suspected of being multiple due to its peculiar 
near infrared spectrum. A compilation of the observed 
properties of all of the sources is provided in Table 1. 

2.2. Imaging and Data Reduction 

Table 2 provides a log of our HST observations. Each 
target was imaged over one orbit in the three filters 
F090M (except SDSS 1254-0122), FllOW and F170M 
using the highest-resolution camera NICl (pixel scale 
0'.'043, field of view ll"xll"). 2MASS 0518-2828 was 
also observed with the F145M and F160W filters. The 
FllOW and F170M filters sample the peak spectral flux 
of T dwarfs (around 1.2 /xm) and the 1.6 /zm CH4 
band, respectively, as illustrated in Figure 1. As the 
near infrared CII4 ba nds are primary classifi cation di- 
agnostics for T dwarfs iBurgasser et alJl2'006aD . FllOW- 
F170M color can provide a rough estimate of spectral 
type (cf. § 3.2) as well as a discriminant for bona-fide, 
low-temperature brown dwarf companions. The F090M 
filter samples the red wing of the pressure-broadened 
0.77 nm K I doublet (Burrows. Marlev. fc S harDll2000t 
lAllard e t al. 2003; Bur rows & Volobuvcv 2003), and pro- 
vides an additional discriminant against background 
sources. 

All data were acquired in MULTIACCUM mode. Mul- 
tiple exposures in the FllOW and F170M filters were ob- 
tained in a spiral dither pattern with steps of 1'.'3 (~ 30 
NICl pixels). Total integration times in these two filters 
ranged over 791-912 s and 1519-1600 s, respectively, for 
the majority of our sample. Exceptions include 2MASS 
0348-6022, which was observed for a longer period due 
to its location in the HST continuous viewing zone; and 
2MASS 0518-2828, where shorter exposures were taken 
to allow observations in five filters over one orbit. Short 
(48-88 s) , single F090M exposures were obtained for the 
GO-9833 targets. Multiple F090M, F145M and F160W 
exposures were obtained for 2MASS 0518-2828 using the 
same dither pattern as the FllOW and F170M observa- 
tions. 

Several of our targets were not well centered on the 

We use abbreviated notation for sources in our observed sam- 
ple throughout the text; e.g., 2MASS hhmmiddmm, where the 
suffix is the J2000 sexigesimal Right Ascension (hours and min- 
utes) and dechnation (degrees and minutes). Full source names 
and coordinates are provided in Table 1. 
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NICl caraera due largely to their uncertain or unknown 
proper motions at the time of the observations. The most 
extreme case is that of 2MASS 0727+1710. At the time 
of the HST observations the proper motion of 2MASS 
0727+1710 had not been measured, so the telescope was 
pointed at the 1997.83 epoch position as measured by 
2MASS. Unfortunately, this source has one o f the largest 
proper motions in our sample, 1'.'297±0'.'005 ijVrba et al.l 
[2004), and the resulting 8'.'3 motion between the 2MASS 
and HST imaging epochs was sufficient to move 2MASS 
0727+1710 out of the NICl field of view. Observations of 
SDSS 0151+1244 were also offset due to source motion, 
and the object was imaged in the corner of the NICl 
camera's field of view (1'.'2 from the closest edge of the 
array), limiting the area sampled for companions. The 
remaining sources were detected sufficiently close to the 
center of the NICl array (>2'.'5 from the array edge) 
to provide adequate sampling of separations within the 
resolution of the original discovery surveys (~ l'.'5-2" for 
SDSS and 2MASS). 

Images were reduced by standard pipeline processing 
(CALNICA, Bushouse et al. 1997) using updated calibra- 
tion images and photometric keywords as of August 2004. 
CALNICA reduction includes analog-to-digital correc- 
tion, subtraction of bias and dark current frames, linear- 
ity correction, correction for readout artifacts (the "bars" 
anomaly), division by an appropriate flat field image, 
photometric calibration, cosmic ray identification, and 
combination of MULTIACCUM frames into a single cal- 
ibrated image. Post-CALNICA processing was limited 
to the cleaning of cosmic rays and persistent bad pixels 
by nearest-neighbor interpolation, and the mosaicking of 
the FllOW and F170M dithered image sets (and all five 
filter sets of 2MASS 0518-2828) using the CALNICB rou- 
tine. 

2.3. Resolved Sources 

Subsections (2'.'5x2'.'5) of the reduced F090M, FllOW 
and F170M mosaic images for each source are shown in 
Figure [3 North/east orientations are indicated by ar- 
rows, and the images are scaled logarithmically to high- 
light low flux features. Note the clearly resolved PSFs in 
most of the FllOW and F170M images, resulting in sig- 
nificant structure outside of the core of the PSF including 
first-order Airy rings and diffraction spikes at wider sep- 
arations. Three sources immediately stand out as ob- 
vious doub les. The previously re ported binary SDSS 
0423-0414 l|Burgasser et al.ll2005bl) shows two overlap- 
ping PSFs roughly oriented along a NNE/SSW axis, 
with the northern component appearing to be slightly 
fainter in the FllOW and F170M bands. SDSS 1021- 
0304 also appears to be a close double aligned along a 
ENE/WSW axis, with the western component appearing 
to be fainter at both F090M and F170M, but not FllOW. 
2MASS 1553+1532 is a well-resolved pair aligned along 
a N/S axis, with the southern component appearing 
to be slightly fainter at F090M (it is marginally de- 
tected in this band), FllOW and F170M. This source 
had previously been reported as acandidate binary by 
iBurgasser. Kirkpatrick Sz Brownl l)2002j) . 

In addition to these three sources, PSFs of 2MASS 
0518-2828 and 2MASS 0926+5847 are shghtly elongated 
in the HST images and therefore also appear to be dou- 
ble. These sources are shown in more detail in Figure |2| 



which displays contour plots of the central 0'.'9x0'.'9 re- 
gions of the F090M and FllOW images around 2MASS 
0518-2828 and SDSS 0926+5847, respectively, and equiv- 
alent data for the unresolved source 2MASS 1503+2525. 
2MASS 0518-2828 is shghtly elongated along a N/S axis 
(PSF full width at half maximum of 1.86 pixels, as com- 
pared to 1.57 pixels for 2MASS 1503+2525), with the 
shape of its southern extension indicating a fainter com- 
ponent. AU three F090M images obtained for 2MASS 
0518-2828 show the same elongation in the same orien- 
tation, lending confidence to its reliability. The elonga- 
tion is less obvious in the FllOW and F145M images 
of this source, and marginally detected in the F160W 
and F170M images, presumably because it is obscured by 
the broader PSF at these wavelengths. SDSS 0926+5847 
is clearly elongated along a NW/SE axis (PSF FWHM 
of 2.67 pixels at FllOW, versus 2.06 pixels for 2MASS 
1503+2525), and appears to be more symmetric, suggest- 
ing near-equal brightness components. Again, the same 
elongation is seen in each of the FllOW and F170M expo- 
sures (the source is only marginally detected at F090M). 
We therefore conclude that both systems are resolved 
doubles. PSF fits for all of the doubles are presented in 
§ 4. The remaining 17 targets appear to be single at the 
resolution of the NICl camera. 

3. NICMOS PHOTOMETRY 
3.1. Measurements 

Aperture photometry for all of the sources in our sam- 
ple were measured from the individual calibrated images 
using the IRAF^^ PHOT routine. Various aperture radii 
ranging from 2-20 pixels (0'.'086-0'.'86) about the source 
flux peak were examined, with a common background 
annulus of 20-30 pixels (0'.'86-l'.'3). Integrated source 
count rates were converted to photometric magnitudes 
on the Arizona Vega system {Mvega = 0.02) using the 
photometric keyword parameter PHOTFNU and Vega 
fluxes of 2157.3, 1784.9 a nd 946.1 Jy at F090M , FllOW 
and F170M, respectively ijSchultz et al J 120051) . Individ- 
ual magnitudes from dithered exposure frames were av- 
eraged to derive a single photometric measurement for 
each source. F145M and F160W magnitudes for 2MASS 
0518-2828 were similarly measured using Vega fluxes of 
1197.1 and 1042.6 Jy, respectively. 

Aperture corrections in the FllOW and F170M fil- 
ters were determined from observations of the three 
highest signal-to-noise (S/N) single sources in our sam- 
ple, 2MASS 0348+6022, SDSS 1254-0122 and 2MASS 
1503+2525. Comparison of integrated flux profiles as a 
function of aperture size between these sources demon- 
strates excellent agreement, with deviations of less than 
0.01 mag for apertures wider than 4 pixels. For F090M 
observations, we adopted aperture corrections measured 
from observations of 2MASS 0518-2828, even though this 
source is slightly resolved (all other sources have insuffi- 
cient S/N in this filter). Table 3 lists the aperture cor- 
rections obtained for each filter, corrected from an 11.5 
pixel (0'.'49) ref erence aperture to an infinite aperture us- 
ing values from iSchultz 

IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Photometric measurements are hsted in Table 4. For 
unresolved sources and the marginally resolved pairs 
SDSS 0926+5847 and 2MASS 0518-2828, we report 5- 
pixel (0'.'22) aperture photometry corrected to an infi- 
nite aperture using the values in Table 3. For the re- 
solved doubles SDSS 0423-0414 and SDSS 1021-0304, we 
report 15-pixel (0'.'65) aperture photometry encompass- 
ing both components with no aperture correction. This 
aperture size was chosen as it includes >90% of the hght 
in all three filters and minimizes photometric noise. For 
the well-resolved double 2MASS 1553+1532, we report 
corrected 3-pixel (0'.'13) aperture photometry for each 
component separately. Uncertainties include contribu- 
tions from the scatter of individual measurements (typ- 
ically 1-2%) and in the aperture corrections (<1% for 
FllOW and F170M, 5% for F090M), as well as 5% ab- 
solute calibrat ion uncertainties and 1% zeropoint drift 
ijSchultz et al.lf2005il . The 5% calibration uncertainties, 
which dominate the error budgets for FllOW and F170M 
magnitudes, are highly correlated and reduce to 3% for 
NICMOS colors; e.g., F110W-F170M. These values are 
reported separately in Table 4. For 2MASS 0518-2828, 
we also measured F145M = 15.86+0.05 and F160W = 
15.20+0.05. 

Formal limiting magnitudes for each source field and 
filter were determined by PSF simulation. Scaled PSFs 
of 2MASS 1503+2525 were added onto blank regions of 
the individual F090M, FllOW and F170M exposures and 
checked for visual detection. Reliable detections were 
possible for peak flux scalings of 7, 5 and 5 times the 
background noise at F090M, FllOW and F170M, respec- 
tively. These limits are given in Table 2. There is a slight 
correlation of these limiting magnitudes with telescope 
pointing angle with respect to the Moon, likely the result 
of increased background emission. Several of our sources 
were detected in the F090M exposures at magnitudes be- 
low the formal detection limits, and have appropriately 
poorer S/N. 

3.2. T Dwarf Colors 

The original motivation for the filter set employed 
in this study was to provide adequate color discrimi- 
nation of bona-fide companions from coincident back- 
ground sources, and to determine photometric classifica- 
tions. Figure 01 compares F110W-F170M colors to spec- 
tral type for sources in our sample and unresolved late- 
type L dwarfs obse rved in the 775 T" NICMOS program of 
iReid et aD l)2006aj) . Spectral types are based on optical 
spectroscopy for the L dwarfs and near infrared spec- 
troscopy for the T dwarfs. We also show synthetic colors 
measured from low resolution near infrared s pectra of 
late-type L and T dwarfs from Burg asser et a l. (2006a). 
Earher than type Tl, F110W-F170M color is relatively 
constant at ^1.8 mag, albeit with significant dispersion 
(±0.3 mag) that is larger than the photometric uncer- 
tainties. For subclasses Tl and later, there is a tight 
correlation between spectral type and color. A linear fit 
to the photometric data for unresolved sources yields 

SpT = 7.26- 3.U{FnOW ~ F170M) (1) 

(where SpT(Tl) = 1, SpT(T5) = 5, etc.), with an RMS 
scatter of 0.4 subclasses. Thus, F110W-F170M is a re- 
liable proxy for spectral type in the T dwarf regime. 
F090M-F110W color also provides a gross discriminant 



of T spectral type, as shown in Figure ]E\ These col- 
ors redden from 1.25 to 2.55 over spectral types Tl to 
T8, due largely to increased absorption by K I. How- 
ever, photometric uncertainties are much larger for the 
F090M data, reducing its utility. We therefore focus on 
the F110W-F170M colors for our analysis. 

Despite the apparent utility of HST colors to distin- 
guish and classify T dwarfs, the vast majority of pho- 
tometric data for these objects are from ground-based 
studies based principally on the J (1.2 /im), H {1.6 /im) 
and K (2.0 ^m) telluric opacity windows. To put our 
photometry into context with existing data, we com- 
pared FllOW and F170M magnitudes to J and H 
photometry, respectively, on the 2MASS and Maun a 
Kea O bservat ory (MKO ; Simo ns &: Tokunag a (2002); 
iTokun aga. Sim ons fc y^^ca (20()3) ') systems. MKO data 
were collated from iLeggett et aLl (|2002); Knapp et alJ 
([200^; and references therein^^. Figure [S] compares J- 
FllOW and i7-F170M colors to F110W-F170M color for 
T dwarfs in our sample. We also plot MKO/NICMOS 
synthetic colors derived from low resolution near infrared 
spectroscopy (from Burgasser et al. 2006a) as a compar- 
ison. There is a marked difference between 2MASS and 
MKO J-FllOW colors. The former are roughly constant 
(—0.75 mag) for -0.2 < F110W-F170M < 1.7, but there 
is significant scatter (^^0.2 mag) due primarily to the 
large uncertainties associated with faint 2MASS T dwarf 
photometry (typically 0.04-0.10 mag). In contrast, MKO 
J-FllOW colors show a tight correlation with FllOW- 
F170M color (MKO photometric uncertainties are typi- 
cally 0.03-0.05 mag) and a 0.2-0.3 magnitude offset from 
2MASS J-FllOW colors due to differences in the filter 
profiles (cf. Stephens & Leggett 2004). A polynomial fit 
to the MKO/NICMOS photometric data for unresolved 
sources yields the relation 

Jmko - FllOW = -1.15 + 0.025(F110M/^ - F170M) 

+0.035iF110W - FUOMf (2) 

with a scatter of 0.04 mag. 2MASS and MKO i7-F170M 
colors are similar (the result of their equivalent _ff-band 
filter profiles), and both span a wider range than J- 
FllOW colors due to strong 1.6 /Ltm CH4 absorption in 
the later-type T dwarfs. 7J-F170M color is also corre- 
lated with F110W-F170M color, and a fit to MKO pho- 
tometry for unresolved sources yields 

Hmko - i^l70M= -0.75 + 0.89(F110iy - F170M) 

-0. 27 {FllOW - F170M f (3) 

with a scatter of 0.05 mag. Interestingly, the combined 
hght colors of SDSS 1021-0304 fall well off both of these 
hues As discussed in § 5.3, these discrepancies are due to 
the unique composition of the SDSS 1021-0304 system. 
Note that equivalent relations for 2MASS/NICM0S pho- 
tometry are presented in Reid et al. ( 2006gi). 

In Figure □ we plot derived FllOW and F170M bolo- 
metric corrections (BC) for T dwarfs in our sample as 
a function of F110W-F170M color. These were com- 
puted from ground-based MKO photometry and i^-band 
BCs from Golimowski et al. (2004|) as, e.g., BCpiiow = 
BCk + Kmko — FllOW . Again, there is good corre- 
lation between the NICMOS BCs and F110W-F170M 

A compilation of these data is maintained by S. Leggett at 
http : //www. jach.hawaii . edu/$\sim$skl/LTdata.html 
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color, and a polynomial fit to the FllOW photometry 
of unresolved sources yields 

BCfiww = 1-43 + 0.11(F110T4^ - F170Af) 

-0.24(F110T4^-F170M)2 (4) 

with a scatter of 0.07 mag. 

4. PSF FITTING 
4.1. Method 

The properties of the five resolved doubles in our sam- 
ple were determined by fitting individual images to PSF 
models using an algorithm similar to that described in 
iBurgasscr et al. ( 2 003c) . One important modification in 
this study was the us e of inodel PSFs generated by the 
Tiny Tim program^"^ ijKristjllQOSl) . Tiny Tim was specif- 
ically designed to generate PSFs for HST imaging data, 
and includes mirror zonal errors and filter passband ef- 
fects in the model PSF shape. We generated several 
grids of Tiny Tim PSFs for the F090M, FllOW and 
F170M filters appropriate for the post-cryocooler NIC- 
MOS NICl detector, samphng 169 positions across the 
chip in row/column increments of 16 pixels. For each 
calibration image, our PSF fitting routine employed the 
model PSF located closest to the position of the target 
on the chip. We also u sed near infrare d spectral data for 
SDSS 1021-0304 from iBureasser et a l. ( 20064 as input 
to the Tiny Tim program in order to model the appropri- 
ate spectral response across the passband. Model PSFs 
were sampled at 10 times the native pixel resolution of 
NICl for subpixel shifting. 

PSF fits were made for 2'.'5x2'.'5 subsections of each 
image centered on the target sources, and initial guesses 
for the pixel positions and fluxes of the two components 
were made using a simple peak detection algorithm (for 
the primary) and single-PSF subtraction (for the sec- 
ondary). Model images were then generated using two 
Tiny Tim PSFs resampled to the resolution of the data 
and scaled to the estimated peak fluxes of the detected 
sources. Our routine then iteratively searched for the 
optimal solution to the primary position, secondary po- 
sition, primary flux and secondary flux, in that order, by 
computing the residuals between the model image and 
the data. Positional shifts of 0.1 pixels were made by 
shifting the oversampled model PSFs in integer units, 
and then downsampling to the resolution of the data. 
Fluxes were varied in steps of 1% (0.01 mag). The PSF 
fits were done recursively for convergence; i.e., if bet- 
ter solutions to any of the four binary parameters were 
found, the routine retested all of the parameters starting 
with the primary position until no reduction in residuals 
could be made. 

Figure IHl illustrates the quality of these fits for one of 
the F170M images of the SDSS 1021-0304 pair. Shown 
are surface plots of the original data on a logarithmic 
vertical scale (to bring up the background noise), the 
best-fit PSF model, the result of subtracting the primary 
PSF model from the data, and the result of subtract- 
ing the full PSF model from the data. The primary- 
subtracted image shows a well-resolved secondary com- 
ponent, with clear detection of that component's first- 
order Airy ring. The final subtraction is extremely clean 



and average residuals (typically <1% of the peak source 
flux for all fits) are at the level of background noise. 

Final estimates of the flux ratio, separation and posi- 
tion angles of each double were determined as the mean 
of individual determinations from each calibrated image, 
with some vetting of very poor fits caused largely by 
cosmic ray hits close to the target source. Uncertain- 
ties include scatter in the individual fits and systematic 
uncertainties of 0.01 mag and 0'.'004 in flux ratio and sep- 
aration, respectively, as prescribed by the fitting routine. 
These values are listed in Tabled 

For the wider doubles SDSS 0423-0414, SDSS 1021- 
0304 and 2MASS 1553+1532, we estimate that addi- 
tional systematic effects in the fitting process are insignif- 
icant given the well-resolved nature of these sources and 
minimal residuals. Indeed, resolved aperture photome- 
try for the 2MASS 1553-1-1532 pair are consistent with 
the PSF results within measured uncertainties. For the 
closer binaries 2MASS 0518-2828 and SDSS 0926-f5847, 
whose angular separations are less than 2 NICl pixels, 
systematic effects may be more important. To examine 
this possibility, we experimented with PSF fits to test 
data constructed to mimic the measured properties of 
the 2MASS 0518-2828 and SDSS 0926-1-5847 pairs in the 
FllOW and F170M bands. A total of 300 test images 
were constructed for each target/filter simulation using 
random pairings of 2'.'5x2'.'5 subsections of images for 
the brightest unresolved sources in our sample (2MASS 
0348-6022, SDSS 1254-0122 and 2MASS 1503-1-2525, for 
a total of 23 PSF images at FllOW and 13 PSF im- 
ages at F170M). The selected PSF images were shifted 
by subpixel resampling to replicate separations and posi- 
tion angles randomly drawn from Gaussian distributions 
centered at the measured values of the binary under in- 
vestigation, and with distribution widths twice the mea- 
sured uncertainties. The secondaries in each test image 
were scaled to a fiux ratio randomly drawn from a uni- 
form distribution (in magnitude space) spanning mag 
to the measured magnitude difference plus three times 
the measured uncertainty. These test images were then 
run through the same PSF fitting algorithm as described 
above to derive experimental values. Systematic effects 
were ascertained by selecting only those test cases where 
the experimental values agreed with the measured val- 
ues for the binaries, and then computing the mean and 
standard deviation of the associated input parameters. 
For 2MASS 0926-1-5847, these simulations indicate that 
the underlying flux ratios for this system (AFllOW^''*™^ 
= 0.4±0.2, AF170M("™) = 0.4±0.3) are closer to unity, 
as suggested by visual inspection of the images them- 
selves (cf. Figure OJ. Similar systematic flux ratio offsets 
are indicated for 2MASS 0518-2828 (AFllOW^"™) = 
0.8±0.5, AF170M(''™) = 0.9±0.6). The flux ratio offsets 
are largely due to the PSF fitting algorithm attempting 
to fit both components with a single PSF, while the sec- 
ondary PSF fits the largest peak in the residuals. There 
were no indications of systematic offsets in the separa- 
tions or position angles of these systems, however. The 
"systematics-corrected" fiux ratios resulting from these 
simulations are given in Table |S1 
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The angular separations of our sources range from 
0'.'051±0'.'012 (2MASS 0518-2828) to 0'.'349±0'.'005 
(2MASS 1553+1532). Assuming spectrophotomet- 
ric distance estimates for 2MASS 0518-2828, SDSS 
0926+5847 and 2MASS 1553+1532 of 34±6, 38±7 
and 12±2 pc based on their combined hght 2MASS 
J-band magnitudes (corrected for equal-brightness 
componen ts) and the Mj / spectral _Wpe relations of 
[Tinnev. B urgasser fc Ki rkpatrick (200^, projected sepa- 
rations range over 1.8-5.0 AU. These values are consistent 
with the small separati ons typically found for resolved 
brown dwarfs (Burgass er et al.ll2006hD . 

Relative magnitudes and combined light photometry 
were used to determine component FllOW and F170M 
magnitudes and colors; these are also listed in Ta- 
bleEl For SDSS 0423-0414, SDSS 1021-0304 and 2MASS 
1553+1532, the secondary colors are consistent with T 
dwarf spectral types of T2, T5 and T7, respectively. 
Component colors for SDSS 0926+5847 (after correc- 
tion for systematic effects) indicate spectral types of 
T4±0.5, consistent with the composite spectral type of 
T4.5. Component colors for 2MASS 0518-2828 have 
much larger uncertainties, and we can only ascertain 
that they are consistent with spectral types '^T3 and 
earlier. In all cases, we can rule out that the secon- 
daries are hotter background stars, since both compo- 
nent and composite F090M magnitudes would be signif- 
icantly brighter than observed. Furthermore, the likeli- 
hood of a background source lying near any of the target 
sources is very small. In our entire sample, only five 
additional point sources were detected at FllOW in the 
11"^ NICl field of view, with magnitudes of 19.5-22.5. 
Assuming that the background surface density scales as 
IQp.QFiww scahng a.S(P), then the probability of a 
background source with FllOW < 18.0 mag (bracketing 
the estimated magnitudes of the detected secondaries) 
lying within 1" of any target source is 5xl0~®, and can 
be ruled out at the >4cr confidence level. We therefore 
conclude that all five secondaries are physically bound T 
dwarf companions. 

Utilizing our derived FllOW bolometric correc- 
tion/color relation (Eqn. 4), we can determine the rel- 
ative bolometric luminosities of the binary components 
as 

/\Nh,i=A-h,i{B)~Ah,i{A) 

= /\FIIQW + BCfiww{B)-BCfiww{A)IS>) 

These values are given in Table |S1 We also list absolute 
Mboi values for the components of SDSS 0423-0414 and 
SDSS 1021-0304, which have parallax distance measure- 
ments. In all cases, we verify that the secondaries are less 
luminous than the primaries, as expected. Relative ef- 
fective temperatures, Tg/Tyi, were determined from the 
relative bolometric luminosities assuming identical com- 
ponent radii, so that Tb/^a — {Lb / La)'^^^ ■ Again, sec- 
ondary Te//S are less than primary Tg^js for all five 
systems, although in no case are differences more than 
~20%. These ratios are consistent with estimated com- 
ponent Te/ fS based on their spectral types (from FllOW- 
F170M color and spectral decompo sition; see jj 5.3) 
and the Tg/ / /spectral type relation of lGolimowski et al.l 
(|2b04') ■ taking into account the 124 K scatter in the lat- 
ter relation. For SDSS 0423-0414 and SDSS 1021-0304, 



we derived component Te//S from their individual lumi- 
nosities and assumed radii of 0.095+0.010 R0, appro- 
priate for 0.5-5 Gyr brown dwarfs i n the Tpff range of 
late- type L and T dwarfs (jBurrows et al.lll997ll . These 
effective temperatures are consistently 100-200 K lower 
than those based on the iGolimowski et al. (2004) rela- 
tion. Although the deviations are comparable to the 
uncertainties in both this relation and our Te// deter- 
minations, the Golimowski et al. (20*3) results may be 
overestimated in this spectral type regime due to con- 
tamination by these previously unresolved (and hence 
overluminous) binaries. 

Finally, we derived mass ratios for the five binary 
systems, assuming co evality, using the ma ss-luminosity 
power-law relation of iBurrows et a D ll200l . L oc M , 
implying 

g«10^n5^. (6) 

System mass ratios are all 0.7 or greater, similar to most 
currently known brown dwarf pairs. Individual compo- 
nent masses were est i mat ed from the evolutionary mod- 
els of IBurrows et al.l l|l997f) using the component Te//s 
(for 2MASS 0518-2828, SDSS 0926+5847 and 2MASS 
1553+1532) or Ahois (for SDSS 0423-0414 and SDSS 
1021-0304), and assuming an age range of 0.5-5 Gyr, typ- 
ical for local disk dwarfs (Rcid & Hawlev 2000). These 
component masses are roughly consistent with the esti- 
mated mass ratios. Orbital period estimates were derived 
assuming circular orbits and semimajor axes a « 1.26p 
ijFischer fc Marcvlll992(l . These periods range over 10- 
50 yr, with the SDSS 0423-0414, 2MASS 0518-2828 and 
SDSS 0926+5847 systems (P < 20 yr) appearing to be 
the best targets for dynamical mass measurements. 

4.3. Search Limits for Very Faint Companions 

In order to search for even fainter companions to 
unresolved T dwarfs in our sample, we repeated the 
PSF fitting analysis described above on the FllOW 
calibrated images using a single PSF model for the 
primary. The FllOW images were chosen because 
this filter samples the peak of the near infrared spec- 
tral fiux of brown dwarfs down to Tg// « 500 K 
(jBurrows. Sudarskv fc Lunindl2003() . and therefore pro- 
vides the most sensitive probe for low mass companions. 
We examined residual images (subtraction of the PSF 
model from the data) by eye for faint point sources that 
persist in the same location relative to the original pri- 
mary. Out of the entire sample, only one "candidate" 
companion was identified, a faint source (FllOW w 22.0) 
located 0'.'97 northeast of the Tl SDSS 0151+1244. This 
source has a high probability of being spurious, how- 
ever; it does not appear in either the F090M or F170M 
images, and due to the pointing offset of the HST ob- 
servations of SDSS 0151+1244 (see § 2.2) was only de- 
tected in one FllOW exposure. While it is therefore 
likely to be a residual cosmic ray or ghost, it warrants 
follow-up confirmation imaging since its brightness rel- 
ative to SDSS 0151+1244 is consistent with a Te// ~ 
450 K brown dwarf companion. 

Sensitivity limits for putative faint companions were 
quantified by examining the FllOW radial profiles for 
each of the unresolved sources before and after PSF sub- 
traction. Examples of the brightest (2MASS 1503+2525) 
and faintest (SDSS 0207+0000) sources are shown in Fig- 
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ureO PSF subtraction results in residuals that are '^4- 
4.5 mag fainter than the peak source flux in the core and 
first Airy ring {p < 0'.'2), irrespective of the brightness of 
the primary. At larger separations, there is little or no 
improvement in sensitivity beyond the inherent decrease 
in the primary flux, and residuals are largely background 
limited at p > 0'.'4. Faint source sensitivity limits (as- 
suming 3cr detections) for each of the unresolved sources, 
including mass ratio limits assuming AFllOW ^ AAf(,oi 
and using Eqn. 6 are given in Table These detection 
limits can be characterized as follows: 

• No detections for p < 0'.'04, 

• AFllOW - 3-3.5 mag (q > 0.3-0.4) for 0'.'04 < p< 
0'.'2, 

• AFllOW varying from ~3-3.5 mag to the back- 
ground Umit (-4-6.5 mag; q > 0.1-0.3) for 0'.'2 < 
p < 0'.'4, and 

• AFllOW background hmited for p > 0'.'4. 

Assuming all of our sources have masses below the hy- 
drogen burning minimum mass (—0.075 Mq; Chabrier et 
al. 2000; Burrows et al. 2001), these observations rule out 
companions down to just above the deuterium burning 
limit (—0.012 M©; Burrows et al. 2001) for most of our 
targets. 

Finally, we comment on the 5 faint sources (FllOW = 
19.5-22.5) detected at wider separations {p > 4") in the 
NICl images. These can be ruled out as low tempera- 
ture companions based on their magnitudes and FllOW- 
F170M colors, typically 1.2-1.4 mag, inconsistent with 
co-spatial mid- and late-type T dwarf companions but 
typical for background M stars. We conclude that no 
bona-fide companions are present around any of our tar- 
gets with p > 0'.'4 and g > 0.1-0.3. 

5. ANALYSIS OF INDIVIDUAL SOURCES 

5.1. SDSS 0423-0414 

The identification of SDSS 0423-0414 as a binary sys- 
tem was previously reported in Burgasscr ct al. (2005b). 
The parameters reported here supersede those of the 
previous paper, although all measurements are consis- 
tent within the reported uncertainties. This source, 
similar in composition to the recently resolved binary 
2MASS J22521073-1730134 (Rcid et al. 2006b), is an un- 
usual system, as its combined light optical spectrum 
ijCruz et al.l 120031 Kirkpatrick et al. in prep.) exhibits 
both 6563 A Ha emission, an indicator of magnetic ac- 
tivity; and 6708 A Li I absorption, present in brown 
dwarfs with masses below the Li-burning minimum mass 
(-0.065 M^: IRebolo"Martm. fc Maeazzu (1992)). Both 
signatures are rare in very late-type L dwarfs and T 
dwarfs, as the strength and frequency of Ha emis- 
sion plummets across the L dwarf regime (Gizis et ajj 
12000; Kirkpatrick ct al. 2000; Mohantv & Basri 2003 
IWest et al.ll2004f) , while Li I absorption becomes increas- 
ingly difficult to detect against a progressively fainter 

The most interesting source is located 7'.' 63 northeast of SDSS 
1254-0122 at 12''54'"54!14 -01''22'42'.'84. With F110W-F170M = 
2.31it0.14, this source is likely to be a highly reddened background 
star or faint unresolved galaxy. 



continuum s uppressed by pressure- broadened Na I and 
K I lines lIBurrows VolobuvevI The only 

other sources known to exhibi t the same combination 
of fea tures are the L2 Kelu 1 l)Ruiz. Leggett. fc AllardI 
'1997), whi ch has a lso been resolved as a binary sys- 
tem (L iu fc Leggetti .20Q5: Gelino. Kulk arni fc Stephen? 
2006); the L0.5 2MASS J2057540 9-0252302 llCruz et a" 
I2003i|) . which has not (|Reid et al.. .2006aj) : and the LO 
2MASS J11544223-3400390 (Kirkpatrick et al., in prep.), 
which has not yet been imaged at high angular resolu- 
tion. 

This raises the question: from which component or 
components do these spectral features arise? Decom- 
position of the combined light near infrared spectrum 
using the resolved NICMOS photometry indicates that 
this system is com posed of an L6 .5 primary and a T2 
secondary (Burgas ser et al.ll2005bl see § 5.3). Ha emis- 
sion from T dwarfs is rare; only three other T dwarfs have 
been detected in emission (Burgasscr ct al. 2003 ai), one 
of wh ich, 2MASS J12373919+6526148 ( Burgasser et alJ 
Il999|) . is unusually active and is speculated to be a 
very tight (p — 0.1 R©) in teracting binary system 
l|Burgasser et al.l2000all2002a|l . The Ha line flux as mea- 
sured from combined light optical spectral data (Kirk- 
patrick et al. in prep.), flux calibrated to SDSS i' pho- 
tometry (20.22±0.04; Geballe et al. 2002), is 1.7x10"^^ 
ergs cm~^ s~^. Using the component bolometric lumi- 
nosities listed in Table [3 we derive log LHa/Lboi = -5.5 
if the emission arises from the L6.5 primary, and -5.2 if it 
arises from the T2 secondary. Compared to similar-typed 
objects exhibiting Ha emission ~ e.g., DENIS-P J0205.4- 
1159 (L7; log LHa/Lboi < -6.2; Mohanty & Basri 2003) 
and SDSS 1254-0122 (T2; log Lna/Lboi = -5.8; Burgasser 
et al. 2003a) - the emission flux from either component 
is not necessarily extreme, but is nevertheless rare in 
this spectral type regime (cf. Figure 3 in Burgasser et al. 
2002a). 

Turning to the 6708 A Li I line, iLiu fc Leggett (|2005fl 
have pointed out that the detection of this feature in a 
substellar binary can be used as a powerful constraint of 
systemic age in conjunction with theoretical evolutionary 
models, particularly if the absorption can be attributed 
to one or both components. SDSS 0423-0414 exhibits a 
prominent Li I line, and as luminous flux in this spec- 
tral region is dominated by the earlier-type primary (the 
optical classification of this source is L7.5; Cruz et al. 
2003) it is likely that this component is responsible for 
the absorption. This deduction is supported by the fact 
that atomic Li gas is likely to be depleted in the pho- 
tosphere of the secondary (for which we derive Tg// = 
1260±70 K) as Li is incorporated into LiCl ar id LiOH 
at Tfi ffS < 1500 K and pressures > 1-10 bar (iLodderj 

Assuming then that the observed absorption arises 
from the L6.5 primary, which must then have M < 0.065 
Mq, an upper age limit of 1.7 Gyr can be deduced us- 
ing th e theoretical evolutionary models of lBurrows et alj 
(^13); this is illustated in Figure EH This age is on the 
young side for a field dwarf, and may explain the presence 
of Ha emission in one or both components; observations 
of fuUy-convective lower main sequence stars in the field 
and clusters shows that magnetic activity i s commonly 
enhanced in young stars ijHawlev et al.ll999l:lReid.200a) . 
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It is also consistent with the kinematics of this system, as 
its tangential velocity, Vtan — 24.0±0.7 km s ~^ , is on the 
low end of the T Dwarf Vtan distribution of iVrba et al.l 
(|2004). An estimated minimum age for this system of 
0.5 Gyr can be argued from the absence of low surface 
gravity features in its combined-light optical spectrum 
(e.g., VO absorption and weakened alkali lines; see Kirk- 
patrick 2005). Thus, the components of SDSS 0423-0414 
are among the few brown dwarfs with well-constrained 
ages, distances and bolometric luminosities. 

5.2. 2MASS 0518-2828 

The detection of a companion to 2MASS 0518-2828 
appea rs to confirm the binary hypothesis of iCruz et al.l 
((2004) for this source, put forth to explain its unusual 
near infrared spectrum. 2MASS 0518-2828 exhibits clear 
CH4 absorption at 1.6 /j,m but no CH4 band at 2.2 /im. 
This is in contrast with trends in the standard L/T spec- 
tral sequence, where the 2.2 /im band is seen to develop in 
the latest- type L dwarfs first, follow ed by 1.6 /i m absorp- 
tion at the start o f the T sequence ((Geballe et a l. 2002) . 
iCruz et ahl l(2004f l found that the combination of L6 and 
T4 spectra, with the latter scaled to be 20% brighter at 
1.27 /im adequately matches the observed spectral energy 
distribution for 2MASS 0518-2828. The small separation 
of this source, and the corresponding poor determina- 
tion of its relative photometry, makes it impossible to 
verify the conjectured spectral types of the components 
based on these ifST observations (although our photom- 
etry are consistent with these types). However, the fact 
that 2MASS 0518-2828 is resolved into two components 
makes this scenario likely. 

5.3. SDSS 1021-0304 

Like SDSS 0423-0414 and 2MASS 0518-2828, SDSS 
1021-0304 appears to be a binary straddling the L/T 
transition, with component F110W-F170M colors indi- 
cating spectral types of <T2 and T5. However, this 
source is particularly interesting as the two components 
have nearly equal magnitudes at FllOW, while the sec- 
ondary is a full magnitude fainter at F170M. Because 
H2O and CH4 absorption bands encompassed by the 
FllOW filter bandpass (cf. Fig. 1) become stronger with 
later spectral types, the equivalent magnitudes of the two 
components suggests that the secondary has a brighter 
peak flux density, hkc 2MASS 0518-2828. 

To explore this possibility, we performed a spectral de- 
composition of the combined light near infrared spec- 
tra of this source and SDS S 0423-0414 using a me thod 
similar to that described in lBurgasser et aTl l|2005b() and 
Held et al. (2006b). In brief, our technique involves 
the combination of various pairings of standard spec- 
tra (sources with well-defined classifications) after scaling 
them to the relative fluxes of the binary system under in- 
vestigation. The hybrid spectra were then compared to 
the combined light spectrum of each (unresolved) binary 
to determine the best match. We performed our analysis 
on low resolution (A/AA ^ 150) near infrared data ob- 
tained with the SpeX spectrograph (Ravncr et al. 200^ 
mounted on the 3.0m NASA Infrared Telescope Facility. 
Details on the acquisition, reduction an d characteristics 
of thes e data a re described in detail in Burgasse r et al.l 
Il2004al I2n06ar . The comparison basis set was com- 
posed of equivalent spectra of late L and T dwarf stan- 



dards from iKirkpatrick et alJ (|1999D ; ICruz et all (|2003(1 ; 
and [Burgasser et al. {2000^; specifically: 2MASS 
J08354256-0819237 (L5), 2MASS J04390101-2353083 
(L6.5), DENIS-P JO205.4-II5915 (Delfosse et al. 199"! 
L7), 2MASS J16322911+1904407 (L8), SDSS 0151+1244 
(Tl), SDSS 1254-0122 (T2), 2MASS 1209-1004 (T3), 
2MASS 2254-f3123 (T4), 2MASS 1503-^2525 (T5) and 
SDSS 1624+0029 (T6). Pairings of the standard spec- 
tra were scaled to the observed FllOW flux ratios, then 
added together and normalized. The quality of agree- 
ment between the resulting hybrid spectra and those of 
the binaries was quantitatively determined by compari- 
so n of the H^O , CH4 and K/J spectral indices defined 
in lBurgasser et al. ( 2006^, as well as the relative F170M 
flux ratios between the spectral components. 

The best matches for SDSS 0423-0414 and SDSS 1021- 
0304 are shown in F igure ITTI For the former, we confirm 
previous results by iBurgasser et al.l (|2005bfl . finding a 
best fit to a combination of the L6.5 2MASS 0439-2353 
and the T2 SDSS 1254-0122. A hybrid spectrum of the 
Tl SDSS 0151+1244 and the T5 2MASS 1503+2525 pro- 
vides the best match for SDSS 1021-0304. Note that the 
derived spectral types of the secondary components agree 
with photometric classifications based on F110W-F170M 
colors. In both cases the hybrid spectra show remark- 
able agreement with the binary spectra, both in terms 
of band strengths and the overall spectral energy distri- 
bution, across the full 0.8-2.5 /im band. The relative 
F170M magnitudes as measured from the scaled compo- 
nent spectra are in rough agreement with ifST photome- 
try, differing by at most ~0.1 mag, an offset attributable 
to the low resolution and calibration uncertainties in the 
spectral data. 

Examining the relative fluxes of the best-fit component 
spectra for SDSS 1021-0304 in more detail, a remarkable 
fact is revealed. The emergent flux density of the T5 
secondary of this system is 31% brighter than that of 
the Tl primary at the peak of the spectral energy dis- 
tribution (1.27 /im), and 24% brighter at the 1.05 /im 
flux peak. This is despite the fact that the secondary is 
11% cooler and 37% less luminous overall. That these 
spectral peaks differ significantly between the two com- 
ponents while FllOW magnitudes are roughly equivalent 
can be explained by the redistribution of flux within the 
FllOW bandpass. The increased brightening at 1.05 and 
1.27 /im in the secondary component is offset by deeper 
H2O and CH4 bands at 1.1 and 1.35 /im. At the bottom 
of these molecular features, and at shorter and longer 
wavelengths, the primary component is brighter. The 
components of SDSS 0423-0414 do not exhibit this same 
brightness inversion, although the peak flux densities are 
rather similar (differing by less than 25% at 1.27 /im) 
given the large differences in spectral type. The bright- 
ening of the secondary component of SDSS 1021-0304 is 
similar to that hypothesized for 2MASS 0518-2828, and 
more recently observe d for the T1.5+T 5.5 binary SDSS 
J153417.05+161546.1 Hhm et a,].l '200fi. hereafter SDSS 
1534+1615). Indeed, the photometric and spectroscopic 

15 DENIS-P J0205. 4-1159 is either a resolved, near-equal 
mass binary (Koerner et al. 1999; Leeeett et al. 2002) or a triple 
IBouv et al.i2005.V and is arguably a poor choice for this analysis. 
However, we did not have an alternate L7 comparison source, and 
DENIS-P J 0205. 4-1159 is currently the L7 optical spectral stan- 
dard on the IKirkpatrick et al.l 119991) scheme. 
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properties of SDSS 1021-0304 and SDSS 1534+1615 are 
quite similar. We discuss the observed J-band bright- 
ening, and its imphcations on the transition between L 
dwarfs and T dwarfs, in § 7. 

5.4. 2MASS 1217-0311 

iBurgasser et al.l ()2003cD reported a possible faint com- 
panion to this object in HST Wide Field Planetary 
Camera 2 imaging, but our NICMOS observations fail 
to reveal this source. Assuming that the relative flux 
at FllOW between 2MASS 1217-0311 and the putative 
companion is as bright or brighter than that at F1042M 
(Ac = 1.02 fiui), it would have been easily detected at 
the separation (0'.'21) and flux ratio (2.4 mag) previously 
observed. It is likely that the prior detection was an un- 
fortunate combination of cosmic ray hits localized near 
the target source. 

5.5. 2MASS 1553+1532 

2MASS 1553-1-1532 is the latest-type binary in our 
sample, and the best resolved. The F110W-F170M 
colors are consistent with very similar spectral types 
(T6.5 and T7), effective temperatures (within 7%) 
and masses {q = 0.90±0.02). With a separation of 
0'.'349±0'.'005 the 2MASS 1553-^1532 pair is resolvable 
by ground-based imaging under the best seeing condi- 
tions, and was previously reported as a p ossible binary by 
IBurgasser. Kirkoatrick fc Brownl 1)20021) based on imag- 
ing observations with the Keck 10m Near Infrared Cam- 
era iMatthews & Soifcr 1994, hereafter NIRC) on 2000 
July 22 (UT). We have revisited these data to determine 
whether the two components share common proper mo- 
tion, and to search for orbital motion. 

Conditions during the Keck observations in 2000 were 
particularly excellent, with clear skies and seeing of 0'.'3 
(full width at half maximum) at Kg during the obser- 
vations. 2MASS 1553-1-1532 was observed in this fil- 
ter, with 10 dithered exposures of 20s each obtained se- 
quentially. Immediately following these observations, 10 
dithered 20s exposures of the unresolved T dwarf 2MASS 
2254-1-3123 were also obtained. Despite the large angular 
offset between the sources, we used these observations for 
PSF calibration as no other sufficiently bright and unre- 
solved sources were detected in the 2MASS 15534-1532 
fields. Raw images for both datasets were pairwise sub- 
tracted to eliminate first order background emission, and 
checked for linearity. No further reduction of the data 
(e.g., fiat fielding) was done as relative photometry was 
not a priority. Figure [HI shows the PSFs of 2MASS 
2254+3123 and 2MASS 1553+1532. The latter is clearly 
extended along a NNE/SSW axis, but the underlying 
components overlap substantially. We extracted astro- 
metric information using a PSF fitting algorithm sim- 
ilar to that described above, but in this case compar- 
ing 3" x3" subsections of each pair- wise subtracted frame 
of 2MASS 1553+1532 to all 10 observations of 2MASS 
2254+3123, for a total of 100 separate fits. Imposing the 
condition that average residuals be less than 3% of the 
peak source fiux, the 42 best fits gave a mean separation 
p = 0'.'30±0'.'03 (assuming a camera pixel scale of 0'.'153) 
and position angle 6 — 199±7°. 

The difference in epoch between the HST and Keck 
images is 3.126 yr. A preliminary proper motion of this 



object as measured by the USNO infrared parallax pro- 
gram l|Vrba et al.ll2004D of ~0'.'4 yr"! (F. Vrba 2006, pri- 
vate communication) implies a total motion of the sys- 
tem of ~1'.'3 over this period. Yet the change in the 
relative separation between the two components is only 
0'.'05±0'.'03. These observations therefore confirm com- 
mon proper motion for this pair, which are almost cer- 
tainly gravitationally bound. 

As for orbital motion, both the separation and position 
angle of this system have changed only slightly between 
the Keck and HST observations. While only marginally 
significant {Atheta = 9±7°), Figure IT!^ illustrates that 
these slight changes are perceptible. However, the small 
position angle change is significantly less than that ex- 
pected (^^26°) for its estimated 45 yr orbital period. This 
suggests that the true orbital period may be much longer, 
possibly due to smaller component masses in a younger 
system, or a particularly eccentric orbit; or that the or- 
bital inclination is quite different from a face-on projec- 
tion (consistent with the slight change in the separation 
of the two components) . Further high-resolution imaging 
may constrain these possibilities in a reasonably short 
time period (~5 yr), but mapping of the full orbit is 
clearly a long-term prospect. 

6. AN UPDATED ASSESSMENT OF BROWN DWARF 
MULTIPLICITY 

6.1. The Binary Fraction 
The fraction of resolved binaries in our sample is ef"^ 
= 5/22 = 23lg^%, where the u ncertainties take into ac - 
count the size of the sample (I Burgasser et all f2003c^) . 
This is similar to resolved binary fractions measured 



field dwarfs llReid et alJ 12001 


I2006al IBouv et al 




IClose et al.l 120031 IGizis et al.l 


20031 ISiegler et al.l 


20051). 



However, as all of these samples are largely magnitude- 
limited, this fraction is biased in favor of unresolved, 
near-equal mass ratio binaries. We can estimate the un- 
derlying (i.e., vo lume-limited) binary fra ction, £(,, as (cf. 
Eqns. 4 and 5 in IBurgasser et all l|2003d) 'l 



^obs 



a(l 



^obs \ 



^obs ' 



where 



^ /,^(l + g^-64)3/2/(g)dg 

f{q)dq 



(7) 



(8) 



is the fractional increase in volume sampled for binaries 



with flux ratio fs/ f. 



.,2.64 



(Eqn. 6) and r nass ra - 



tio ratio distribution f{q). In ' Burgasser et all 
both flat and delta- function forms of /(g) were consid- 
ered; here, we explicitly calculate a — 2.50tQ og using 
a power-law mass ratio distribution as described below. 
This yields e?, = 11^3%, a value consistent with previous 
determinations of bias-corrected VLM binar y fractions 
(.Burgasser et al,. .2Q03c; Close et al, 2003: Siegler et alJ 
12005ft and volume-limited estimates ()Bouv et al J 12003 
iReid et al.ll2006al -15% and 12^^%, respectively). 

It is important to stress that this fraction applies only 
to those binary systems that are resolvable by direct 
imaging. For our study, this limits the phase space sam- 
pled to p > 1 AU and q > 0.35, or p > 5 AU and q > 
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0.2. The equivalent ph ase space of F-G and M dwa rf bi- 
naries in the studies of iDuauennov fc Mavod l)199lD and 
IReid fc Cxiziij iITqqI yield binary fractions of 39% (com- 
bining both mass ratio and period distributions) and 
24^4% (for My > 9), respectively. Hence, in equiva- 
lent phase spaces the binary fraction of brown dwarfs in 
our sample is less than that of more massive stars. 

But does this mean that the overall binary fraction 
of brown dwarfs is less? Maxtcd & Jeffries (2005) have 
proposed that a substantial fraction (50-67%) of VLM 
binaries may be hiding in more closely separated sys- 
tems {p < 2.6 AU) and can only be resolved as spectro- 
scopic binaries. This projected separation corresponds 
to an angular separation of <0'.'13 at the average dis- 
tance of sources in our sample ('^20 pc), only 3 NICl 
pixels. Indeed, over 25% of resolved brown dwarf bi- 
naries identified to date have angular separations below 
this limit, with the majority close to the resolution limits 
of HST. This supports the possibility that a significant 
number of more closely-separated and/or more distant 
systems remain unresolved. Bayesia n statistical analysi s 
of high resolution imaging studies bv IReid et a l] ll2006ar . 
which takes into account the possible presence of unre- 
solved systems, indicates an overall VLM binary fraction 
of 24%, twice that of the resolved fraction. This still 
places the binary fraction of VLM dwarfs at 1 /3 that of 
solar-type stars, and 2/3 that of M dwarfs, consistent 
with a decreasing binary fraction toward later spectral 
types and lower masses. 

6.2. The Separation Distribution 

The projected separation distribution of 30 brown 
dwarf binaries resolved by high resolution imaging sur- 
veys to date are shown in Figure HH These include sys- 
tems listed in Burgasser et al. (2006b)^^ that have esti- 
mated primary masses below 0.075 Mq. In accordance 
with previous studies, we find that this distribution peaks 
at very close separations, '^4 AU with a broad peak 
spanning 2-8 AU. This is significantly lower than the 
30 AU pe ak of the F-G and M binary separation dis- 
tributions ijDuauennov &: Mavor l'1991: Fisch er fc Marcvl 
|l992). Indeed, no brown dwarf field binaries have been 
identified with separations > 15 AU. However, two wider 
brown dwarf binaries systems have been recently iden- 
tified in young cluster/ associations: 2MASS Jl 101 1926- 
7732383 ^Luhm an"20041 . a 240 AU binary in the -2 Myr 
Chameleon I associat ion; and 2MASS J1207334-393254 
ijChauvin et all 12004 |l005), a 40 AU, very low mass 
(Mtot ~ 0.03 Mq) br own dwarf b inary in the ~8 Myr TW 
Hydrae association l)Gizisll20'o3) . Both systems are very 
young, and it remains unclear as to whether their con- 
figurations are stable long-term ( Mugrauer fc Neuhiiuseil 
|2Q05). The widest VLM fiel d binary so far ide ntified, DE- 
NIS J055146.0-443412.2 (Bi Ueres et al.ll2005l 220 AU), is 
likely composed of two VLM stars. Hence, the wide sep- 

1® Note that this analysis assumes a symmetric Guassian sepa- 
ration distribution; the fraction may be higher (lower) if there is an 
excess (deficiency) of short period systems (cf. Maxted & Jeffries 
2005). 

^"^ T his sample incorpor a tes binaries identified in |Martin_et_aI 

e;'Leinert et al." f200l');'Reid et al."f200l','2006a); 'Potter et al 
; Bouv ct al. 12003); Burgasser et al. 1 2003c); Gizis ot al 
: iMcCaughrean et al . (2004); .Burgasser et al.. l.2005aF : 
ILiu fc Leggeta 120051) : ILiu et al. (200^; and this study. 



ar ation brown dwarf b inary desert originally suggested 
bv iMartin et alJ l)2000D — not to be confused with the 
brown dwa rf companion desert around solar-type stars 
ijMarcv fc B utler 2000) — remains a distinct character- 
istic of brown dwarf binaries in the field. 

At closer separations, imaging surveys are limited by 
angular resolution. Hence, the true peak of the brown 
dwarf separation distribution may be lower than that in- 
ferred from Figure ITTI Nevertheless, it is interesting to 
note that the separations of binaries identified in this sur- 
vey — and indeed all T dwarf binaries identified to date 
— are < 5 AU, such that most T dwarf binaries have sep- 
arations below the peak of the brown dwarf distribution. 
The apparent compactness of T dwarf binaries as com- 
pared to warmer M- and L-type brown dwarf systems is 
consistent with a maxi mum binar y separation that scales 
with total system mass Reid et a l. 2001; Burg asser et alJ 
lc1loseet~aLll2003r since cooler brown dwarfs have 
lower masses than warmer ones at a given age. However, 
the smaller separations of T dwarf binaries may also be 
due to selection effects. T dwarfs are intrinsically fainter 
and typically found at closer distances to the Sun than 
M and L dwarfs in magnitude-limited surveys. Hence, T 
dwarf binaries can generally be observed at higher linear 
resolution. A statistically robust volume-limited sample 
of M-, L- and T-type brown dwarfs would provide an 
adequate check for mass dependency in the separation 
distribution of substellar objects. 

6.3. The Mass Ratio Distribution 

The mass ratio distribution of brown dwarf binaries is 
shown in Figure IT^ This distribution is clearly peaked 
at q ~ 1, with 50±9% of all known systems having near- 
equal mass components. Again, because the majority of 
these systems were originally selected from magnitude- 
limited surveys, there is an inherent bias in the discov- 
ery of equal-mass systems that scales approximately as 
(1 -I- qr2-64-j3/2_ bias-corrected distribution, also shown 
in Figure IT^ nevertheless shows that near-equal mass ra- 
tio systems are predominant. This result is robust even 
when sensitivity limits are taken into account. The bias- 
corrected frequency of binaries drops by a factor of 7.8 
from q = I to q — 0.5, even though most imaging pro- 
grams are complete for companions down to or below this 
limit. A fit to the bias-corrected distribution for q > 0.5 
to a power-law, f{q) oc g^, yields 7 = 4.2±1.0, slightly 
flatter but neverthele ss consistent with a Bayesian anal- 
ysis of VLM binaries ijReid et alJl2006aj) . 

In summary, our sample supports prior results on 
brown dwarf multiplicity, namely: 

• The resolved binary fraction of brown dwarfs is 
lower than that of stars, Cb ~ 11% for p > 3 AU 
and q > 0.3; 

• The separation distribution of resolved brown 
dwarfs peaks around 4 AU; the true peak may lie at 
lower separations due to resolution limits of imag- 
ing programs; 

• The maximum separations of field brown dwarf bi- 
naries appears to decrease for later spectral types, 
consistent with a mass-dependent trend; and 
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• Most brown dwarf pairs have near-equal mass ra- 
tios, with a bias-corrected distribution of f{q) oc 
g,(4.2±i.o) infjicated by current data. 

These characteristics of brown dwarf field binaries pro- 
vide key empirical constraints for the theoretical mod- 
eling of brown dwarf formation and dynam ical evolu- 



tion, iss u es that are discu s sed in detail i n.Burgasser et al 
|20061)l) : iLuhman et al] l)2006|) : and IWhitworth et al 

7. BINARIES AND THE L/T TRANSITION 

7.1. J-band Brightening - Intrinsic to the L/T 
Transition 

Three of the binaries in our sample ~ 2MASS 0518- 
2828, SDSS 0423-0414 and SDSS 1021-0304 - are com- 
posed of brown dwarfs that span the transition between 
L dwarfs and T dwarfs. This spectral type range has 
been the focus of both observational and theoretical stud- 
ies as it encompasses dramatic changes in the atmo- 
spheric properties (e.g., photospheric condensate dust 
depletion) and spectral energy distributions (e.g., the 
onset of CH4 absorption) of cool brown dwarfs. This 
transition also exhibits several unusual traits, includ- 
ing an apparent brightening of absolute J-band magni- 
tudes from late-type L to mid-type T dwarfs jPahu et al, 
' 20021 iTinnev. Biirgasser fc Kirknatric k 2003 rVrba et all 



20041) . This so-called "J-band bump" has been at- 
tributed to dynamic atmosph eric processes, such as c on- 
densate cloud fragmentation fBurgasser et al."20021 ^. a 
sudden increase in sedimentation cfhcicncy ( Kna pp et al.l 
I2OOI or a glob al collapse of the condensate cloud 
layer ijTsuiil 120051) . However. iTsuii fc Nakaiimal l|200l 
have also argued that age and/or surface gravity effects 
amongst disparate field sources may be responsible. 

The component fluxes of the SDSS 1021-0304 binary 
demonstrate that the last hypothesis can be largely ruled 
out for this system. Under the reasonable assumption of 
coevality, these brown dwarfs have similar ages, masses 
and (presumably) radii, implying nearly identical surface 
gravities. Yet the T5 secondary of this system is clearly 
brighter than the Tl primary at 1.05 and 1.27 /im. Simi- 
lar trends suggested in the 2MASS 0518-2828 system and 
observed in SDSS 1534-1-1615 demonstrate that SDSS 
1021-0304 is not a unique case. Hence, a brightening 
of surface fluxes at these wavelengths appears to be an 
intrinsic feature of the L/T transition. 

7.2. A "Bump" or a "Plateau"? 

In th eir analysis of the SDSS 1534-1-1615 binary, 
iLhi et al. ( 2006) proposed that the J-band bump may 
be artificially enhanced by a significant contribution 
of binaries amongst mid-type T dwarfs (such "crypto- 
binarity" has also been suggested by Burrows, Su- 
darsky & Hubeny 2006). To examine this hypothesis 
in detail. Figure ^| compares absolute MKO J- and 



dwarf systems with measured parallaxes \ 


Dahn et al. 


2002 


; iTinnev. Burgasser & KirkDatrickll2003l 


Vrba et al. 


2004 


, precision >20%), companions to nearby Hippar- 



cos s t ars llBe gklii^fc Zuckcrman 1988; Nakaiima et^D 



1^95: Burgass er et all I2000bt iKirkpatrick et alJ |2I, 
McCaughrean et al.ll2004D and resolved absolute magni- 



tudes for SDSS 0423-0414, SDSS 1021-0304 and the Tl 



+ T6 binary Epsilon Indi B l|McCaughrean et al.ll200^ . 
For consistency, spec tral types for L d\yarfs ar e based on 
optical data and the IKirkpatrick et al"! l)1999l) classifica- 
tion scheme, while tho se for T dwarfs a re based on near 
infrared data and the 'Burgasse r et al.l (j2006a) scheme. 
MKO J magnitudes for the SDSS 0423-0414 and SDSS 
1021-0304 components are based on their F110W-F170M 
colors and Eqn. 2; JiT-band component photometry is de- 
rived from synthetic colors measured from the compo- 
nent spectral templates. 

The J-band bump is seen clearly in these data in 
the T1-T5 spectral type range, and is well traced 
by the absolute m agnitude /spec tral t ype relation of 
ITinnev. Burgasser fc Kirkpatrickj l|2003j) . Yet one of 
these data points is the unresolved SDSS 1021-0304 
systems, and its individual component fluxes are only 
slightly brighter (J = 14.33 and 14.29) than the latest- 
type L dwarfs (J ~ 14.7). The same holds true for the 
T2 secondary of SDSS 0423-0414 (J = 14.38) and the 
Tl primary of Epsilon Indi B ( J = 14.30). Indeed, all of 
the resolved components spanning types Tl to T5, have 
nearly identical absolute J-band magnitudes. 

Does this mean that the J-band bump is largely an ar- 
tifact of multiplicity? Possibly, but only if the T2 SDSS 
1254-0122, the T3.5 SDSS 1750-^1759 (both unresolved 
in this study) and the T 4.5 2MASS J0 5591914-1404488 
llBurgasser et al.l l2000d unresolved in IBurgasser et alJ 
l)2003d) ') are all closely-separated multiples. This is 
not out of the realm of possibility, for as discussed in 
§6.2 the separations of brown dwarf binaries likely ex- 
tend below imaging resolution limits. One or all of 
these systems may also have been imaged at an un- 
f ort unat e orbital angle, as was t he case initially for 
Kelu 1 llMartm. Brandner fc Basri||1999l iLhi fc Leggettl 
120051 iGehno. Kulkarni fc Stephen &^'2006'). Furthermore, 
the fact that SDSS 1254-0122 and 2MASS 0559-1404 are 
^0.5 and '^0.8 mag brighter at J-band than the T2 and 
T5 secondaries of SDSS 0423-0414 and SDSS 1021-0304, 
respectively, suggests that the former are unresolved mul- 
tiples. If the primaries of these hypothetical systems are 
constrained to have Mj « 14.3, then the secondaries of 
SDSS 1254-0122 and SDSS 1750+1759 would have Mj « 
15.3-15.5 and be ~T6 dwarfs. 2MASS 0559-1404 would 
be required to have near equal-magnitudes components 
or be a higher multiple system. High resolution radial ve- 
locity monitoring observations are needed to test these 
possibilities. 

It is important to point out that the absolute J-band 
magnitudes of the early-type T dwarf resolved binary 
components examined here are still ~0.4 mag brighter 
than the latest-type L dwarfs (but ~0.5 mag fainter at 
JC-band). Hence, some broad-band brightening may still 
be present across the L/T transition. Furthermore, even 
if the J-band "bump" is a more modest "plateau" , there 
remains a significant ('^30%) brightening at 1.05 and 1.27 
/im due to flux redistribution within the J-band spectral 
region, a feature not yet reproduced self-consistently by 
current atmosphere models. 

7.3. The Origin of J-band Brightening 

That the observed brightening is concentrated in the 
1.05 and 1.27 /im flux peaks is an important clue to 
its origin. The photospheric atomic and molecular gas 
opacities of low-temperature brown dwarfs show distinct 
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minima at these wavelengths, shaped by strong H2O and 
CH4 bands and bracketed by pressure-broadened K I at 
shorter wavelengths and collision-induced H2 absorption 
at longer wavelengths. Condensate opacities, for grain 
sizes (^40-80 /im) computed in a self-consistent manner 
ifAckernian fc Marlev 200L1 . are roughly constant across 
the near infrared band. In the L dwarf regime, conden- 
sates are a dominant source of opacity at the J-, H- 
and i^-band spectral peaks. However, theoretical at- 
mosphere models incorporating condensate clouds indi- 
cate that the photospheric opacity from these species are 
dominant only in the 1.05 and 1.27 /zm flux peaks for 
Teff < 1300-1500 K; i.e., at the L/T transition (cf. Fig- 
ure 16 in Burrows, Sudarsky & Hubeny 2006). If these 
condensates are suddenly removed, the total opacity at 
these wavelengths decreases, allowing brighter emission 
from deeper and hotter layers. 

This is the underlying thesis for dynamical atmo spheric 
explanations for the J-band brightening ( Bureasse r et al.1 
I2002bl: iKnapD et'al] l2004() . However, one must also 
consider whether higher gas opacities at longer wave- 
lengths, with the increased photospheric abundances of 
H2O and CH4 molecule s below T^jf ^ 1300-18 00 K 
ijBurrows fc Shard IT999t iLodders fc E^glevl I200I and 
stronger H2 absorption, might lead to a redistribution 
of flux into the 1.05 and 1.27 /im flux peaks. Current 
cloud models that assume constant sedimentation effi- 
ciencv l|M arley et alj|2002 D or p article size distributions 
ijBurrows. Sudarskv fc Hubenvl I2OO6I do not show this 
to be the case. Hence, a dynamic mechanism for clear- 
ing out photospheric condensate dust may still be neces- 
sary to explain the evolution of brown dwarf atmospheres 
across the L/T transition. 

7.4. The Frequency of L/T Binaries 

Is there evidence from the distribution of binary 
frequency as a function of spectral type that bi- 
naries play a special role through the L/T transi- 
tion? To address this, we have compiled results 
from high res o lution imaging of L and T dwarfs by 
IKoerner et al." |I19 99D: iMartm et alJ lll999H: I Reid et al.' 
jlpOl, 200 6a): C lo se et all (120031): IBouv et a l. (2003); 
iB^rgasser et'alJ ()2003c^) : iGizis et al.l ()2003j) : and this 



study. We did not incl ude individual binary discover- 



ies made serendipitously 


Goto et alJl2002l IPotter et all 


120021 iMcCauehrean et al. 


120041 iBureasser et al."2005a() 



or those identified a s part of as yet unpublished sur- 
veys pu et al. 2006": St umpf. Bradner fc Henninsll20d^ 



in order to make a fair assessment of the observed bi- 
nary fraction. Care was taken to identify duplicate 
sources in each of the imaging studies, and classifications 
wer e verified through published optical (L dwarfs, on 
the iKirkpatrick et al. (1999) sc heme) and nea r infrared 
data (T dwarfs, on the Burgasser et al.l ()2006a|) scheme). 
Only those resolved pairs that had a high probability of 
companionship, based on either common proper motion 
confirmation, resolved spectroscopy and/or photometric 
colors, or very low probability of coincidence with an un- 
related background source, were considered as bona-fidc 
binaries. The complete sample includes 129 L dwarfs and 

We also include the recently identified Kelu 1 binary 
ILiu fc Leggetai200a ILtelino. Kulkarni fc tjtephensll200ai . a target 
of prior searches. 



34 T dwarfs, of which 33 are binary. 

Figure [T7I plots the observed binary fraction of these 
sources as a function of spectral type, binned by indi- 
vidual subclasses and into subclass groups of L0-L2 (62 
sources), L2.5-L4.5 (28 sources), L5-L6.5 (27 sources), 
L7-L9.5 (12 sources), T0-T3.5 (7 sources), T4-T5.5 (12 
sources) and T6-T8 (15 sources). Note that these frac- 
tions have not been corrected for selection bias (resulting 
in an overestimate from equal-brightness systems) or sen- 
sitivity/resolution limits (resulting in an underestimate 
by missing closely separated or low mass ratio systems). 
This sample may also be subject to more subtle biases, 
such as the smaller typical distances of later-type, intrin- 
sically fainter brown dwarfs, resulting in greater linear 
resolution for these objects (although this effect may be 
offset by the apparent decrease in separations for lower- 
mass brown dwarfs). It nevertheless serves to illustrate 
possible trends. 

There is clearly significant structure in the binary frac- 
tion distribution for individual subclasses, although this 
could be attributed to small number statistics. By bin- 
ning the subclasses (reducing statistical uncertainties), 
a remarkable result emerges. For most of the sample, 
binary fractions are consistent with the overall fraction, 
ef" = 20±4%. Yet the L7-L9.5 and T0-T3.5 subclass 
groups - the L/T transition objects — have fractions that 
are twice as high, 42^![q% combining all 19 systems in 
this spectral type range. This deviation is significant at 
the 98% confidence level compared to the sample mean. 

Why would the observed binary fraction of L/T transi- 
tion objects be so high? We posit the following scenario. 
Analysis of the SDSS 0423-0414 and SPSS 1021-0304 
compo n ents, and prior r esults f ro m IKirkpatrick et alJ 
(12000): iB^reasser et all (I2002d): IDahn et alJ (|2002|) : 
IVrba et all (1 2004*) and 'Goh mowski et all (|20041 . aU in- 
dicate that the L/T transition spans a relatively nar- 
row range of effective temperatures, ATeff ~ 200- 
300 K. However, the cooling rate of brown dwarfs is 
larg ely insensitive to cha nges in the photospheric opac- 
ity (jChabrier et al.ir2000 ^. such as the removal of con- 
densates or emergence of CH4 absorption. Brown dwarfs 
must therefore progress through the L/T transition rel- 
atively rapidly, implying fewer such sources per spectral 
subtype for a given field sample. On the other hand, 
the analysis of §5.3, and similar results by 'Cruz et alJ 
(12004): Bu rgasser ct al. (2005b); Rcid ct al.'( 2006b); and 
ILiu et allf2006j) . all demonstrate that early- type T dwarf 
spectral features can be reproduced from the combined 
light of a late-type L and mid-type T dwarf binary. It is 
therefore possible that such hybrid binaries, if unrecog- 
nized, could significantly contaminate a spectral sample 
of early- type T dwarfs. 

To illustrate this point, consider the following example. 
Assuming that L5-L8 dwarfs have T^ff « 1700-1300 K 
and L8-T5 dwarfs have T^ff ~ 1300-1100 K (Table 5 
andiGolimowski ct al. (2004)), the mass function simula- 
tions of lBurgasseil ()2004D predict a relative space density 
of Nl/t/Nl « 0.9 between these two groups, largely 
independent of the shape of the underlying mass func- 
tion. However, because mid- to late-type L dwarfs are 
roughly twice as bright as L/T transition objects, the 
relative number observed in a magnitude- limited sample 
(the best approximation for current imaging samples) is 
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(Nl/t/Nl)"'"' « 0.3. Now consider that all brown dwarfs 
in a magnitude-limited sample have a resolvable binary 
fraction of ~25%, and that 20% of all late-type L dwarf 
binaries have T dwarf secondaries (this is roughly consis- 
tent with the mass ratio distribution of Figurc ll5|l . These 
binaries would exhibit a combined light spectrum simi- 
lar to a late- type L/early-type T dwarf, and would be 
identified as such in an unresolved spectroscopic sample. 
Hence, the observed binary fraction among late-type L 
dwarfs in this scenario would be ^^20%, while the fraction 
of L/T transition binaries would be 



obs 

'l/t 



0.2x0.25Nl + 0.25Nl/t _ 0.05 + 0.075 
0.2x0.25Nl + Nl/t ^ 



0.05 + 0.3 



36%; 

i.e., nearly twice the underlying binary fraction. This 
numerical example serves to illustrate that the binary 
hypothesis provides both a qualitative and quantitative 
explanation for the peak in the binary fraction of L/T 
transition objects. More complete modeling of this effect 
will be presented in a forthcoming publication. 

We therefore conclude that multiplicity does play an 
important role in the L/T transition, contaminating sam- 
ples of "true" transition objects and leading to a greater 
J-band brightening than that inferred for resolved sys- 
tems. These binaries also provide a detailed and intrigu- 
ing picture of this still poorly-understood transition, and 
a list of all currently known L/T binaries is given in Ta- 
ble [7| Further study of these source will provide im- 
proved understanding of the physical mechanisms gov- 
erning this transition, including the depletion of photo- 
spheric condensates, the emergence of CH4 gas and the 
possible role of atmospheric dynamics in brown dwarf 
spectral evolution. 

8. SUMMARY 

We have identified 5 binaries in of a sample of 22 T 
dwarfs imaged with HST NICMOS. Of these, three are 
well-resolved, permitting determination of their compo- 
nent spectral types, relative bolometric luminosities and 
TgffS, and systemic mass ratios. The identification of 
2MASS 0518-2828 as a closely-separated binary confirms 
previous suspicions of multiplicity based on this object's 
unusual near infrared spectrum. The bias-corrected re- 
solved binary fraction of this sample (e;, = 1X^3%), the 
near-unity mass ratios of the components of these sys- 
tems {q > 0.7) and their small projected separations 
(p < 5 AU) are all consistent with previously identified 
trends amongst VLM dwarfs, indicating that they are 
salient properties of brown dwarf field binaries. 

Three of the binaries in our sample, SDSS 0423-0414, 
2MASS 0518-2828 and SDSS 1021-0304, are composed 
of sources spanning the L to T transition, and spec- 
tral decomposition analysis of SDSS 1021-0304 reveals 
that its T5 secondary is 25-30% brighter at 1.05 and 
1.27 fim than its Tl primary, despite being 35% less lu- 
minous overall. The properties of these sources, as well 
as the recently discovered binary SDSS 1534+1615, in- 
dicate that the J-band brightening previously observed 
amongst late-type L and mid-type field T dwarfs is an in- 
trinsic feature of the L/T transition and not the result of 
age, surface g ravity or metallicity effects. In suppo rt of 
the results of iBurrows. Sudarskv fc Hubenvl l)2006D and 
iLiu et all 1)2006(1 . we find that the J-band bump may be 



more of a J-band "plateau" , with T1-T5 dwarfs having 
Mj « 14.3, enhanced by the presence of unresolved bi- 
naries in this spectral type range. Indeed, we find that 
the frequency of L/T transition binaries is twice as high 
as those of all other L and T dwarfs, a statistically sig- 
nificant deviation that can be explained if spectroscopic 
samples of L/T transition objects are significantly con- 
taminated by binaries composed of earlier-type and later- 
type components. Taken together, the properties of L/T 
binary systems provide the most conclusive evidence to 
date that the L/T transition occurs relatively rapidly, 
driven by the removal of photospheric condensates that 
is likely to be facilitated by dynamic atmospheric pro- 
cesses. Further parallax and multiplicity measurements 
will better constrain the flux evolution and relative num- 
bers of L/T transition objects, important constraints for 
understanding the physical mechanism of photospheric 
condensate depletion and the atmospheric evolution of 
brown dwarfs as they cool below Tg// « 1500 K. 
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22'' 54™ 


18?92 


+31°23'49'.'8 


1998.48 


15.26+0.05 


15.02+0.08 


14.90+0.15 






2MASS J23391025+1352284 


T5 


23''39™ 


10?25 


+13°52'28'.'4 


2000.91 


16.24+0.11 


15.82+0.15 


16.15+0.31 




0.83i 



Refere nces. — 



Dahn et, all 



) Gcballo ct al. 12002): (2) 
(7) .Burgasscr ct al.. 1.2003h 



BurKasscr ct al. (2002c); (3) Burgasscr. McElwain & Kirkpatrick (2003); 
; (8) .Strauss ct al.. I.199al : (9) .Vrba ct al.. I.2004h C 10) .Tinncv. Burgassci 



(4) Cruz ct al. (2004); (5) Leg 
& Kirkpatricta ^2003.) : fll) .Tin] 



^Near infrared spectral types from lBurgasser et alj i2006aD . 



Coordinates from the 2MASS All Sky Point Source Catalog 
'^Discovery reference in boldface type, followed by references for additional photometric and astromctric data. 
<*MKO JHK from lLeggett et al.l I2002H or lKnapp eTal] 1200411 . 



TABLE 2 

Log of HST Observations for Programs GO-9833 and GO-10247. 



F090M FllOW F170M 



Object 


UT Date 


t 




t 




t 




Roll Angle' 






(s) 


(mag) 




(mag) 


(s) 


(mag) 


(°) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


SDSS 0151+1244 


2003 Nov 3 


56 


19.0 


864 


22.9 


1519 


Zl , 


,Z 


212.0 


SDSS 0207+0000 


2004 Feb 3 


48 


18.4 


864 


22.1 


1519 


on 


o 
.y 


208.1 


2MASS 0243-2453 


2004 Feb 15 


56 


19.5 


864 


23.6 


1519 


91 


.0 


220.6 


2MASS 0348-6022 


2004 May 26 


88 


20.2 


1599 


22.7 


3071 


91 


A 


321.1 


2MASS 0415-0935 


2004 Feb 20 


56 


19.0 


864 


22.6 


1519 


91 

Zl , 


1 

, i 


214.1 


SDSS 0423-0414 


2004 Jul 22 


72 


19.3 


816 


22.6 


1519 


91 

Z_L , 


Q 
.O 


15.3 


2MASS 0516-0445 


2004 Jul 25 


48 


18.6 


864 


22.7 


1519 


21, 


.1 


15.9 


2MASS 0518-2828'^ 


2004 Sep 7 


960 


21.0 


64 


19.9 


416 


19, 


.9 


38.9 


2MASS 0727+1710 


2004 Mar 27 


56 


19.5 


864 


23.6 


1519 


21, 


.1 


232.7 


2MASS 0755+2212 


2004 Mar 29 


56 


18.6 


864 


22.4 


1519 


20, 


.9 


236.9 


SDSS 0837-0000 


2004 Feb 18 


48 


18.6 


864 


23.3 


1519 


21, 


.7 


196.5 


SDSS 0926+5847 


2004 Feb 5 


48 


19.1 


960 


22.7 


1600 


21, 


.2 


319.5 


SDSS 1021-0304 


2004 May 22 


48 


18.9 


864 


22.9 


1519 


21, 


.0 


245.9 


SDSS 1110+0116 


2004 May 22 


48 


18.9 


864 


23.3 


1519 


21, 


.5 


247.1 


2MASS 1217-0311 


2004 Apr 26 


48 


18.7 


864 


22.7 


1519 


21, 


.3 


255.7 


SDSS 1254-0122 


2004 Feb 13 






896 


22.7 


1519 


21, 


.1 


64.9 


2MASS 1503+2525 


2003 Sep 9 


80 


18.9 


791 


21.9 


1519 


20, 


.7 


216.0 


2MASS 1553+1532 


2003 Sep 7 


56 


19.5 


864 


23.5 


1519 


21, 


.7 


226.9 


SDSS 1624+0029 


2003 Sop 9 


48 


18.8 


864 


22.7 


1519 


20, 


.9 


230.8 


SDSS 1750+1759 


2003 Sop 12 


56 


18.6 


864 


22.1 


1519 


20, 


.9 


231.8 


2MASS 2228-4310 


2004 May 26 


72 


19.2 


912 


23.0 


1519 


21, 


.3 


37.3 


2MASS 2254+3123 


2003 Sop 15 


64 


19.4 


864 


22.8 


1519 


21, 


.2 


299.7 


2MASS 2339+1352 


2003 Nov 1 


56 


18.8 


864 


22.6 


1519 


20, 


.9 


218.5 



^Limiting magnitude estimated for a 5cr flux peak detection. 

'^Telescope roll angle, East from North. 

''Also observed at F145M (320 s) and F160W (160 s). 



TABLE 3 

NICMOS NICl Aper,tur,e Corrections. 



Aperture 



Radius 


FOgOM'' 


FllOW 


F170M 


(pixels) 


(mag) 




(mag) 


(mag) 


2.0 


-0.71+0, 


,05 


-0.79+0.02 


-1.096+0.009 


2.5 


-0.60+0, 


.04 


-0.695+0.017 


-0.869+0.012 


3.0 


-0.47+0, 


.04 


-0.637+0.014 


-0.757+0.008 


3.5 


-0.35+0, 


.04 


-0.557+0.011 


-0.714+0.008 


4.0 


-0.26+0, 


.04 


-0.451+0.008 


-0.694+0.007 


4.5 


-0.22+0, 


.04 


-0.349+0.007 


-0.658+0.007 


5.0*^ 


-0.20+0, 


.05 


-0.275+0.005 


-0.593+0.006 


5.5 


-0.18+0, 


.05 


-0.232+0.005 


-0.503+0.006 


6.0 


-0.17+0, 


.04 


-0.211+0.004 


-0.410+0.006 


6.5 


-0.16+0, 


.04 


-0.199+0.003 


-0.333+0.006 


7.0 


-0.15+0, 


.03 


-0.190+0.003 


-0.279+0.005 


7.5 


-0.14+0, 


,03 


-0.181+0.002 


-0.247+0.004 


8.0 


-0.14+0, 


.03 


-0.1736+0.0019 


-0.232+0.004 


8.5 


-0.14+0, 


,04 


-0.1681+0.0015 


-0.225+0.004 


9.0 


-0.13+0, 


,04 


-0.1638+0.0012 


-0.219+0.003 


9.5 


-0.12+0, 


.04 


-0.1597+0.0010 


-0.212+0.002 


10.0 


-0.12+0, 


.04 


-0.1549+0.0007 


-0.2039+0.0017 


10.5 


-0.11+0, 


,02 


-0.1498+0.0006 


-0.1969+0.0013 


11.0 


-0.115+0, 


,012 


-0.1448+0.0006 


-0.1919+0.0009 



Note. — Values include corrections from a reference aperture 
(11.5 pixels) to an infi nitely sized ape rture of -0.1136, -0.1393 and 
-0.1888 mag as given in lSchultz et alj 1200511 . 

^Measured from observations of 2MASS 0518-2828. 
^Adopted aperture radius for single sources. 



TABLE 4 
NICMOS Photometry. 



Object 


SpT 


F090M 


FllOW 


F170M 


F110W-F170M 




(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



bubb U151+lz44 


Tl 


i5.5o±(J.U9 


17.2d±U.U5 


15.o4±U.U5 


i.7z±U.U4 


CT^cc non'7 i f\rif\f\ 
bJJob U2U/+UUUU 


T4.5 


1 n n_Ln o 


i / .o4±U.Uo 


Ib.o ( ±U.Uo 


U.y ±U.U4 


OA/T A GC nO/1Q OA^Q 

ziVlAoo U^4o — Z4iOO 


1 D 




iD.ZoitU.UO 


1 c; '7 A -Un Oft 


U.4yztU.U4 


zMAob Uo4o— dU2z 


T7 


-| o cr 1 j_n no 


iD.li±U.UD 


lb. iDzbUAJo 


-U.UOzbU.I.Jo 


zMAoo U4io— UyoO 


T8 


iy.U4±U.iU 


lb. 4 r ±U.Uo 


Ib.D ±U.Ub 


n on _L n n /1 
-U.2U±U.U4 


oJJoo U4zo— U414 


1 L) 


iD.OOltU.iU 


io.zo±u.uo 


io.bz±U.Uo 


l.bD=tU.U4 


zivi/\oo uoiD — yj^^o 


1 O.O 


1 Q Q«_i_n fiQ 
io-yoniu. uy 


iD.ODniU-UO 


ID. l-rtlCU-UO 


fl e:nif\ r\o 
U. 0ZniU-Uv5 


2MASS 0518-2828'' 


Tip 


18.55+0.08 


16.69+0.07 


14.92+0.05 


1.77+0.06 


2MASS 0755+2212 


T5 


18.89+0.09 


16.58+0.05 


16.08+0.05 


0.50+0.03 


SDSS 0837-0000 


Tl 


>18.6 


17.91+0.05 


16.21+0.05 


1.70+0.04 


SDSS 0926+5847 


T4.5 


18.66+0.09 


16.57+0.05 


15.64+0.05 


0.93+0.03 


SDSS 1021-0304 


T3 


19.34+0.15'' 


17.09+0.05'' 


15.83+0.05'' 


1.26+0.04" 


SDSS 1110+0116 


T5.5 


19.15+0.12 


17.27+0.05 


16.49+0.05 


0.78+0.04 


2MASS 1217-0311 


T7.5 


18.91+0.09 


16.73+0.05 


16.77+0.05 


-0.04+0.04 


SDSS 1254-0122 


T2 




15.66+0.05 


14.10+0.05 


1.56+0.03 


2MASS 1503+2525 


T5 


16.84+0.08 


14.65+0.05 


14.18+0.05 


0.47+0.03 


2MASS 1553+1532A 


T7 


19.16+0.16 


17.16+0.05 


17.01+0.05 


0.15+0.04 


2MASS 1553+1532B 


T7: 


>19.5 


17.46+0.05 


17.40+0.05 


0.06+0.04 


SDSS 1624+0029 


T6 


18.46+0.09 


16.28+0.05 


15.96+0.05 


0.32+0.03 


SDSS 1750+1759 


T3.5 


18.80+0.09 


17.18+0.05 


16.04+0.05 


1.14+0.04 


2MASS 2228-4310 


T6 


18.77+0.09 


16.36+0.05 


16.05+0.05 


0.31+0.03 


2MASS 2254+3123 


T4 


18.12+0.10 


16.05+0.05 


15.09+0.05 


0.96+0.04 


2MASS 2339+1352 


T5 


18.86+0.09 


16.91+0.05 


16.27+0.05 


0.64+0.03 



Note. — Aperture photometry derived using an 5-pixel aperture and aperture corrections from 
Table 2, unless otherwise noted. Values are given on the Arizona Vega system (Mvega — 0.02). 
Uncertainties include RMS scatter in count rates between individual exposures, aperture correc- 
tion uncertainties, 5% flux calibration (3% in F110W-F170M color) and 1% sensitivity variation 
(zeropoint drift). 

^Photometry measured using a 15-pixel aperture and no aperture correction to incorporate both 
components. 

^Additional photometry for 2MASS 0518-2828: F145M = 15.86±0.05, F160W ^ 15.20±0.05. 



TABLE 5 
Binary Properties. 



Parameter 


SDSS 0423-0414 


2MASS 0518-2828 


SDSS 0926-1-5847 


SDSS 1021-0304 


2MASS 1553-1-1532 


UT Date 


2004 Jul 22 


2004 Sen 7 


2004 Feb 5 


2004 May 22 


2003 Sep 7 


2000 Jul 22 


InstrumGiitr 


HST NICMOS 


HST NICMOS 


HST NICMOS 


HST NICMOS 


HST NICMOS 


Keck NIRC 


(") 


0'.'164±0'.'005 


0'.'051±0''012 


0"070±0"006 


0'.'172±0'.' 005 


0"349±0"005 


0"30±0'.'02 


(AU) 


2.49±0.07 


1.8±0.5 


2.6±0.5 


5.0±0.7 


4.2±0.7 


3.6±0.7 


e (°) 


19?2±0?8 


189°±8° 


132?9±1?9 


244?6±0?8 


189?9±0?4 


199°±7° 


AF090M 


0.88±0.03 


1.6±0.4 




0.21±0.10 






AFllOW 


0.526±0.015 


0.8±0.5'' 


0.4±0.2=' 


0.06±0.04 


0.31±0.04 




AF170M 


0.820±0.013 


0.9±0.6'' 


0.4±0.3'' 


1.030±0.019 


0.461±0.016 




FllOW (A) 


15.80±0.05 


17.11±0.18 


17.14±0.10 


17.81±0.05 


17.15±0.05 




(B) 


16.33±0.05 


17.9±0.3 


17.54±0.13 


17.87±0.05 


17.47±0.05 




F170M (A) 


14.04±0.05 


15.31±0.19 


16.21±0.13 


16.19±0.05 


16.98±0.05 




(B) 


14.86±0.05 


16.2±0.4 


16.61±0.18 


17.22±0.05 


17.44±0.05 




F110W-F170M (A) 


1.76±0.07 


1.8±0.3 


0.93±0.16 


1.63±0.07 


0.17±0.07 




(B) 


1.47±0.07 


1.7±0.5 


0.9±0.2 


0.66±0.07 


0.03±0.07 




Est. SpT 


L6.5+T2'' 


L6:-|-T4:'' 


T4:-|-T4; 


T1-I-T5'' 


T6.5 


+T7 




0.72±0.13 


0.9±0.6 


0.4±0.2 


0.49±0.13 


0.31±0.12 


Mboi (A) 


15.77±0.16 






16.5±0.7 






(B) 


16.50±0.16 






17.0±0.7 






T2/T1 


0.847±0.013 


O.SliO.ll 


0.91±0.05 


0.894±0.015 


0.931±0.014 


Te//<= (K) (A) 


1490±100 


~1600 


~1330 


1260±210 


~980 


(B) 


1250±80 


~1330 


~1330 


1130±190 


~890 


q 


0.78±0.02 


0.74±0.15 


0.87±0.07 


0.84±0.02 


0.90±0.02 


Est. Mass<i (Mq) (A) 


0.039-0.062 


0.042-0.077 


0.029-0.073 


0.025-0.076 


0.019- 


0.065 


(B) 


0.029-0.051 


0.031-0.074 


0.029-0.073 


0.021-0.074 


0.016- 


0.061 


Est. Period" (yr) 


~19 


~10 


-18 


-50 


-45 



^Estimated flux ratios including systematic effects; see §4.1. 
^Based on spectral decomposition; see §5.3. 

'^Based on IGolimowski et all I2004D T^// /spectral type relation, with the exception of SDSS 0423-0413 and SDSS 1021-0304, where T^ffS 
are derived from measured Mj^^is and assumed radii of 0.095±0.10 R© . 

^Based on ages of 0.5-5.0 Gyr (with the exception of SDSS 0423-0414, wh ere 0.5-1.7 Gyr i s assumed), estimated T^/ys (except SDSS 0423-0414 
and SDSS 1021-0304, where A-f^o? is used) and the evolutionary models of lBurrows et all nWH . 

^Assuming semimajor axes a — 1.26p tiFischer Marcviil992iV 



TABLE 6 

Faint Source Detection Limits for Unresolved Targets. 



0''04 <P< 0'.'2 0'.'2 < p < I'.'O 



Object 


Distance^ 


At'llOW 


q 


AFllOW 


q 






(pc) 


(mag) 




(mag) 




(AU) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


SDSS 0151+1244 


21.3+1.4 


3.1 


0.34 


5.3 


0.16 


0.9 


SDSS 0207+0000 


29+8 


2.8 


0.38 


3.8 


0.27 


1.2 


2MASS 0243-2453 


10.6+0.5 


3.1 


0.34 


4.6 


0.20 


0.5 


2MASS 0348-6022 


~7 


3.1 


0.34 


5.4 


0.15 


0.3 


2MASS 0415-0935 


5.75+0.10 


3.2 


0.33 


6.0 


0.12 


0.2 


2MASS 0516-0445 


~19 


3.3 


0.32 


5.9 


0.13 


0.8 


2MASS 0755+2212 


~20 


3.5 


0.29 


5.8 


0.13 


0.9 


SDSS 0837-0000 


29+12 


3.4 


0.31 


3.7 


0.28 


1.2 


SDSS 1110+0116 


~23 


3.0 


0.35 


5.3 


0.16 


1.0 


2MASS 1217-0311 


11.0+0.2 


3.1 


0.34 


5.3 


0.16 


0.5 


SDSS 1254-0122 


13.7+0.4 


3.5 


0.29 


6.5 


0.10 


0.6 


2MASS 1503+2525 


-9 


3.1 


0.34 


6.5 


0.10 


0.4 


SDSS 1624+0029 


10.9+0.2 


3.1 


0.34 


6.0 


0.12 


0.5 


SDSS 1750+1759 


28+4 


3.0 


0.35 


4.6 


0.20 


1.2 


2MASS 2228-4310 


~13 


3.2 


0.33 


5.7 


0.14 


0.6 


2MASS 2254+3123 


~19 


2.9 


0.36 


5.3 


0.16 


0.8 


2MASS 2339+1352 


~26 


3.5 


0.29 


5.6 


0.14 


1.1 



Note. — Given values are 3c7 detection limits for a single PSF subtraction from 
individual calibration images, separated into the core/first Airy ring region (p < 
0. 2) and the background-dominated region [p > 4). Mass ratio (g) limits assume 
AFllOW ~ AMboi and Eqn. 6. ^ 

^Distance estimates for objects without parallax measurements are based on 
apparent 2MASS J-band magn itudes and the Mj /spectral type relation of 
iTinnev. Burgasser 



TABLE 7 
L/T Transition Binaries. 



Name 


Spectral Types 


Distance'' 


Separation 


Period 


Note 




(A) 


(B) 


(pc) 


(") 


(AU) 


(yr) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


SDSS J042348.57-041403.5 


L6 


T2 


15.2+0.4 


0'.' 164+0'.' 005 


2.49+0.07 


~19 


1,2 


2MASS J05185995-2828372 


L6: 


T4: 


~34 


0'.'051+0'.'012 


1.8+0.5 


~10 


3,4 


2MASS J08503593+1057156 


L6 


T: 


26+2 


0'.'16±0"010 


4.4+0.4 


~43 


5 


2MASS J09201223+3517429 


L6.5 


T: 


~21 


0'.'07±0"010 


1.5+0.5 


-6 


5 


Gliese 337C 


L8 


T: 


20.5+0.4 


0" 53+0'.' 03 


10.9+0.7 


~150 


6,7 


SDSS J102109.69-030420.1 


Tl 


T5 


29+4 


0'.' 172+0'.' 005 


5.0+0.7 


~48 


2,4 


Epsilon Indi B 


Tl 


T6 


3.626+0.013 


0'.' 732+0'.' 002 


2.654+0.012 


~15 


7,8 


SDSS J153417.05+161546.1 


T1.5 


T5.5 


~36 


0'.'110+0'.'005 


3.9+0.6 


~28 


9 


2MASS J17281150+3948593 


L7 


T: 


~23 


0'.'131+0'.'003 


3.0+0.5 


~35 


10 


2MASS J22521073-1730134 


L6 


T2: 


14+3 


0'.' 130+0'.' 002 


1.8+0.4 


-9 


11 



Re ferences. — (1) Bureassor ct al. 12005b); (2) Vrba ct al (2004); (.3) Cruz ct al. ( 200J 
IReid ct aH,200L); ( 6) Burgasser ct al. I 2005a) ; (7) .ESA. a997.1 : (8) JilcCauEhrcan ct al. (20^ 
(lO) IGizis et al.ri200ai ; (ll) TReid et all iSooebl) . 



^Parallax distance measurements from IeSAI il997f) and IVrba et al are given with uncertainties, all others 

are spectrophotomctric distance estimates from the discovery references. 
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Fig. 1. — NICl F090M fyellow V FllOW (blue) and F1 70M (r ed) filter transmission profiles overlaid on the red optical and near infrared 
spectrum of 2MASS 1503+2525 feurgasser et al.ll2003all2004al) . Spectral data are normalized at 1.27 /im. Filter transmission profiles are 
preflight measurements and do not include the NICl detector quantum efficiency or optical element response curves. Key H2O and CH4 
bands present in the spectra of T dwarfs arc indicated. 




Fig. 2.— NICMOS F090M (left), FllOW (center) and F170M (right) images of T dwarfs observed in this study. Subsections of 2'.'5x2'.'5 
in size are shown on a logarithmic scale. Orientations North (arrow) and East (line) are indicated in the center panel. 




Fig. 2. — Continued. 
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Fig. 2. — Continued. 
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Fig. 3.— Contour plots of FllOW images of the central 0'.'9x0'.'9 regions around SDSS 0926+5847 (top left) and 2MASS 1503+2525 
(top right), and F090M images of 2MASS 0518-2828 (bottom left) and 2MASS 1503+2525 (bottom right). Lines indicate isofiuxes of 0.01, 
0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.8, 0.9 and 1.0 times the peak source flux. Orientations (North and East) for each 
image are indicated in the lower right corners. 
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Fig. 4. — F110W-F170M colors versus spectral type for subtypes L6 to T8. L dwarf classifications are based on optical spectra 
(e.g., Kirkpatrick et al. 1999 ) ; T dw arf classifications are based on near infrared spectra (e.g., Burgasser et al. 2006). Photometry from 
this program and lReid etai] 12006al) are shown as solid circles with error bars; multiple sources are encircled. Spectrophotometric colors 
measured from low resolution near infrared data from Burgasser et al. 1 2006a) arc shown as small open circles. A linear fit to the photometric 
data of unresolved objects is indicated by the solid line; ±lcr scatter about this relation arc indicated by the dotted lines. 
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Fig. 5. — F090M-F110W colors versus spectral type for T dwarfs in our sample. Symbols are the same as in Figure HI 
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Fig. 6. — Top: J — FUOW versus FllOW — F170M color for sources in our sample. Solid circles indicate colors on the MKO system, 
solid triangles indicate colors on the 2MASS system. Multiple sources are encircled. Spectropliotometric colors on the MKO system 
measured from low resolution near infrared data arc shown as small open circles. A polynomial fit to the MKO photometry of unresolved 
sources is indicated by the solid line; ±1ct scatter about this relation are indicated by the dotted lines. Bottom: H — F170M color versus 
FUOW — F170M color for sources in our sample. Symbols are the same as in the top panel, and only MKO photometry are shown. A 
polynomial fit to photometry for unresolved sources is indicated by the solid line; ±la scatter about this relation are indicated by the 
dotted lines. The discrepant (multiple) source SDSS 1021-0304 is labeled. 
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Fig. 7.— FllOW and F170M bolomctric corrections (BC) as a function of FUOW - FllOM color. BC values arc based on /C-band 
bolometric corrections from Golimowski ct al. (2004) and MKO j-tT-band ph otometry for sources in our sample. Spectrophotometric BCs 
measured from low resolution near infrared data from lBurgasser et alj I2006af) are shown as small open circles. Polynomial fits to the BCs 
of unresolved sources observed in this program are indicated by the solid lines; ±lcr scatter about these relations are indicated by the 
dotted lines. 




Fig. 8. — Surface plots illustrating PSF fits to F170M imaging data of the SDSS 1021-0304 binary. Shown clockwise from upper left are 
2''5x2''5 subsections of imaging data, the best fitting PSF model, subtraction of the primary PSF model from the data, and subtraction of 
the full PSF model from the data. All four plots are normalized to a common logarithmic scale; average residuals in the final subtraction 
are 0.4% of the peak source flux. 
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Fig. 9. — Sensitivity limits for faint companions around the faintest and brightest unresolved sources, SDSS 0207+0000 (top) and 2MASS 
1503+2525 (bottom). Open circles trace the relative brightness profile with respect to the peak PSF in the original NICMOS image; filled 
circles trace the brightness profile after PSF subtraction. The solid lines trace the radial profile of an oversampled model PSF from Tiny 
Tim. Residuals in the PSF core (0'.'04-0"2) after model subtraction are consistently ~4 mag below the peak primary flux. The far wing 
sensitivity is dominated by background noise. 
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Fig. 10. — Comparison of component bolomctric luminositcs for SDSS 0423-0414 to theoretical evolutionary models. Isomass evolutionary 
tracks arc indicated by lines and labeled (in Solar masses). The 1% Li depletion limit is indicated by the thick line, and effectively lies 
along the 0.065 Mq track for the ages shown. A minimum age for the system of ~0.5 Gyr is based on the absence of notable low surface 
gravity features in the optical and/or near infrared spectrum of this object. A maximum age of 1.7 Gyr is constrained by the detection of 
strong Li I absorption, most likely from the primary component. 
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Fig. 11.— Spectral decomposition of SDSS 0423-0414 (left) and SDSS 1021-0304 (right). The top of each panel shows the best fit spectral 
standards to the primary (red dashed line) and secondary (blue solid line) of each binary, normalized and scaled to the relative FllOW 
magnitudes as measured with NICMOS. The bottom of each panel compares the combined light spectrum of the binary (black solid line) 
to the hybrid spectrum of the two spectral components (purple dashed line), both normalized at 1.27 /im. The agreement between the 
combined light and hybrid spectra is overall quite remarkable. 
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Fig. 12.— Contour plots for one set of NIRC i^^-band images of 2MASS 2254+3123 (left) and 2MASS 1553-1-1532 (right). The displayed 
boxes correspond to a l'.'8xl'.'8 area on the sky, and orientation for both images is indicated by the arrow (pointing north) in the right 
panel. Contour levels corresponding to 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 and 0.9 times the peak source flux are shown. Relative separations 
for the 42 best PSF fits to 2MASS 1553-f 1532 are indicated by the plus symbol (corresponding to the primary) and crosses (corresponding 
to the secondaries), and are consistent with a mean separation of 0'.'30±0'.'02 and position angle of 199±7°. 
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Fig. 13^.— Comparison between NIRC Xs-band (contours) and HST FllOW (greyscale) images of 2MASS 1553+1532. Both images 
are orient^ with north up and e^t to the left, and ang^^ar scale is indicated i^the upper left cornerQThe position angles l^e^een the 
primary and secondary as determined by PSF fitting arc overlaid; the single line and dashed lines correspond to mean and la uncertainty 
as measured from the NIRC images, the double line corresponds to the best fit from the HST images. A slight rotation and expansion of 
the system between the two imaging epochs can be discerned. 
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Fig. 14. — Projected separation distribution (light grey histogram) of 30 brown dwarf binaries identified in high resolution imaging 
studies by Martin et al. 11999); Leinert et al. (2001); Rcid et al. 12001, 2006a); Potter et al. (2002); Bouy et al. (2003); Burgasser et al] 
p003c); Gizis et al. (2003); McCaughrean et al. 12004); Burgasser et al. 1 2005a); Liu & Leggett (2005); Liu et al. (2006); and this study. 
Uncertainties based on counting statistics are indicated by error bars. The distribution exhibits a peak at p 4 AU (logjQ p = 0.6) as 
derived from a Gaussian fit (solid line), although the decline at smaller separations may be the result of resolution limits in the imaging 
studies. Note that T dwarf binaries (dark grey histogram) have smaller projected separations on average as compared to all resolved brown 
dwarf binaries. 
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Fig. 15. — Mass ratio distribution of the 30 brown dwarf binaries from Figure lT4l Uncertainties based on counting statistics are indicated 
by error bars, upper limits are indicated by open circles. Bias corrected values are shown by dashed lines. A power law fit to the bias- 
corrected values, f{q) 

(X<?(4.2±1.0)_ is shown by the thick solid line. Sensitivity limits at close (p > 0''04) and wide (p > 0'.'2) separations 

are indicated. 
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Fig. 16. — Absolute MKO J - (left) and j<'-band (right) magnitude s versus spectral type for field L and T dwarfs with parallax 
measurements fDahn et al." "2002; "I'inncv. Burgasser & Kirkpatrick"200i?; 'Vrba et al."2004') and companions to nearby Hipparcos stars 
JBocklin & Zuckcrman 1988; Nakaiima et al. 1995; Burgasser ct al. 2000b; Kirkpatrick et al. 2001; McCaughrcan et al. 2004). Spectral 
types are based on optical data for the L dwarfs and near infrared data for the T dwarfs. Combined light photometry for known binaries 
are encircled. Primary (red circles) and secondary (blue circles) spectral types and magnitudes for the SDSS 0423-0414, SDSS 1021-0304 
and Epsilon Indi B binaries are indicated. Absolute MKO magnitude/spectral type relations from fTinnev. Burgasser fc Kirkpatrickl 12003?) 
are shown by the solid curves. 
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Fig. 17. — Observed b inary fraction s of L and T dwarfs as a function of spectral t ype. Data were co mpiled from the imaging survey s 
of IKoerner et al J IT99Hh : IMartm et alJ tT99flh : iReid et al.l i200ll Htm^ : IClose eTail EOCi^ -. jBouv et al.l : iBurEfasser et aTT 12003 A : 

lGm^^aLn2003) : and this study. Counting uncertainties are indicated in both panels. The left plot shows binary fractions broken down 
by individual subclasses; upper (zero binaries) and lower limits (all binaries) are indicated by arrows. The right plot groups sources into 
spectral class bins of L0-L2, L2.5-L4.5, L5-L6.5, L7-L9.5, T0-T3.5, T4-T5.5 and T6-T8, with the number of source in each bin labeled. The 
overall observed binary fraction, ejj'" = 20±4%, is indicated by the dashed and dotted lines. 



